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[Plates 22, 23.] 

In a previous paper ■^''' I remarked that there might be reason to suppose that the 
earhest form of a sateUite might not be annular. Whether or not the present inves- 
tigation does actually help us to understand the v^orking of the nebular hypothesis, 
the idea there alluded to was the existence of a dumb-bell shaped figure of equili- 
brium, such as is shown in the figures at the end of this paper. These figures were 
already drawn when a paper by M. Poincare appeared, in which, amongst other 
things, a similar conclusion was arrived at. My paper was accordingly kept back in 
order that an attempt might be made to apply the important principles enounced by 
him to this mode of treatment of the problem. The results of that attempt are, for 
reasons explained below, given in the Appendix. 

The subject of figures of equilibrium of rotating masses of fluid is here considered 
from a point of view so wholly different from that of M. PoiNCAHifi that, notwith- 
standing his priority and the greater completeness of his work, it still appears worth 
while to present this paper. 

The method of treatment here employed is simple of conception ; but it is unfor- 
tunate that, to carry out the idea, a very formidable array of analysis is necessary. 

In the last section a summary will be found of the principal conclusions, in which 
analysis is avoided. 

§ 1. Formulce of Spherical Harmonic Analysis. 

Let there be two sets of rectangular axes, as shown in fig. 1; and let z be 
measured from o to 0, whilst Z is measured from O to o ; let r^ = i:c^ -f- 2/^ + ^^ 
R^ = X'' +Y^ + Z^; and let c = oO. ' 

Then 

. x-^X=0, y^Y=0, z + Z=c (1) 

Let Wi, Wi, denote the solid zonal harmonics of degree i of the coordinates x, y, %, 
and X, Y, Z, respectively. 

Now we shall require to express the solid zonal and certain tesseral harmonics of 

^ * Phil. Trans.,' Part II., 1881, p. 534 
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negative degrees with respect to the origin as solid zonal and tesseral harmonics of 
positive degrees with respect to the origin o, and vice versd; moreover, the results 
will have to be applied to a sphere of radius a with centre o, and to a sphere of 
radius A with centre O. This last clause is introduced in order to explain the intro- 
d action of the symbols a, A^ in this place. 




^ 



o 



The formulas required will be called '^ transference formute," because they are to be 
used in shifting the origin from one point to the other. 

The obvious symmetry of our axes is such that every transference formula from 

O to has its exact counterpart for transference from o to O; thus a second 

symmetrical formula with capital and small letters interchanged will generally be left 

unwritten. When necessary, 6, (f>, will be written for co-latitude and longitude with 

regard to x,y,z; and ©, #, for the same with respect to X, F, Z. 

Then, since 

J?^ = r*^ -|~ c^ — 2rc cos 6^ 

we have the usual expansion in zonal harmonics 

The usual formula for the derivation of the zonal harmonic of negative degree 
i + 1 from l/M is 



il dZJ B'^ B^^^^ 



. (3) 



Hence, on differentiating (2) i times with respect to Z, or, which is the same thing, 
with respect to — z, we have, from (3), 

But 

~ Wk = ^ (^ — 1) . . . (^ — i + 1) -Wu^i = ^"-^ ID^^i . 

Hence 



c 
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In interpreting this formula, it will be noted that, if i is less than k, the term 
vanishes : hence the summation runs from ^=00 to k = i ; it is therefore better to 
write k + i for Jc, and we thus obtain 

This is the first transference formula by which the solid zonal harmonic of degree 
— i — 1 with respect to is expressed as a series of solid harmonics of positive 
degree with respect to 0. The formula (4) includes (2) as the particular case where 
i = 0. The right-hand side of (4) is convergent for r less than a. A similar formula, 
convergent for r greater than a, is easily obtainable, but with this we shall not 
concern ourselves. 

It remains to find the transference formula for certain tesseral harmonics. 

If we put 

p =z ^ {^x -f- y ), . . . . . . . ^ . . (5; 

the general expression for the zonal harmonic is 



il 



Wi-^2,{'—') t /> , . . . • . . . ' (6) 



hfA-m\ 



where the summation extends from A; = to ^ = ^i or i|(* — 1 ). 
From (6) we have 

^ — V (^\h - '^^' L^«-3^ ^^-1 /7\ 



Now, since r^ = 2:^ + 4/9, we have 



7^-p = t{-)H\ 



— ik+l) 4th 

mmmmm 



k + V?.i-2k-2\ ' kf.i-2k\ 



z p 



xjliSO 



, [Ah (^ + 1) - (^ - 27.) (^^ -2k- 1)] , 
-^^ ^ *• {k-\-l).kf.i-2k\ ^ P- • • • (8) 



Subtracting (9) from (8), and simplifying the difference, we have 






^ "^ ^' ''' * -^ '-^%'^2k'^% Jc-¥l-l 



J f^ i ■\. 0. : ' * • » * ' • ' * » » ' » » • tXU/ 
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the last transformation being derived from (7) with i + 2 in place of i, and k 
in place of E 

Differentiate (10) with respect to p, and notice that dr^/dp == 4, and we have 



"?j- 1 



"" dp^ "^^^^'"^^^^ ^ dp^ 



Then^ with ^ + 2 in place of i, 



Now 



mS 


d^'Wi^Q 
dp^ 


— 2 


(2i + 3) 


d 
dp 


Wi 
^2i + l 


1 f ^d'Wi 






1 


W*?^! 4-2 !_ 



<iie'*| + 3 t^%^i + 4 



€^p dp 



2 (2i 4" l)t{;^ 



• « • • • 



r 



Si + 3 ^ 



by (10). 



Differentiating again. 



d^ Wi 



dp^ 7 



3 «,2i + l ^,2i + 5 



f 



do^ 



2(2i + 3)^ 



1 c2' 



-.2i + 5 ^ ^2 



dp' 



^^'^ + 4^7 (11), 



or 



1 d^vJi+c^ # -x^^.-^g 



^,2e+l j^p2 gp2 ^,21-3 



. (U) 



• •••** 



• (12) 



But since p = |(a;^ + y^), it follows that, in operating on a function involving x and y 
only in the form x^ 4* ^^^ 



d 

dm 



d 



rv* 



^ dp ' 



Also 



so that 



Now let us put 



fi^ 



Then 



d^ 



d ^ d 

—~~ 23:: A^ — ' 

dy ^^ dp ' 



2 rf;!> ^ *-^ gp 



d? 6? 



ddf w 



3 5 






(titJO 



1 

2 2 



c^2/? 



-- JL ( /V*'^ ,mmm> <?/* | 



d 
dp 



1 /i t'i 



+ w 



1 (a;2 - /) 
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cf' 



# 



C^aJ^ f^' 



3 



» » » 



g^ = I (^^ « 3^2) 



2\ „. 



fl^ 



(ip*'^ 



and therefore (12) may be written 



%. 









d 
dp 



« « 



• (13) 



4«*!|44»* 



(14) 
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These expressions in (14) are obviously solid tesseral liarmoDics. 

The transference formula required is for ''^' ^ ^ • 
By the formula (4) we have 

operating on both sides by S^, and applying (14), we have 

_g2i-fi^^»^^'+s — -^^f^ i^2lk\~<^* ...... (15) 

Now the general formula (6) for the zonal harmonic shows us that dho^/dp^ is zero 
when X' = 0, 1, 2, 3, and hence S^i(?^ vanishes for the same values of h Thus the 
summation in (15) is from ^ == co to ^ = 4, or, if we write ^ + 2 for h, from co to 2. 
Hence (15) gives 

j>gi+i« ^^^ + 3— ^^f^ ^_2!7^ + 2!\c/ # .... \i^y) 

This is the second transference formula required. 

We observe that the transference of a negative zonal harmonic gives us positive 
zonals, and that tesseral harmonics of the type S^W^-^.^/i^^^'*"^ give us harmonics of the 
type SV;^+2. 

§ 2. The Mutual InJluenGe oftivo Spheres of Fluid without Rotation, 

Imagine two approximately spherical masses of fluid of unit density, with their 
centres at the origins o and respectively, and with mean radii a and A respectively. 

We shall find that each exercises on the other certain forces, one part of which has 
a solid zonal harmonic of the first degree as potential. This part of the force must 
remain essentially unbalanced, in the supposed system, but we shall see hereafter that 
it is balanced by the rotation to be afterw^ards imposed on the system. 

Meanwhile it will be supposed that it is annulled in some way, and we shall content 
ourselves with finding the mutual influence of the spheroids, and the outstanding 
term of the first degree of harmonics. 

Let us assume that the equations, referred to our two origins, of the surfaces of the 
two spheroids, when they mutually perturb one another, are 



T 






/^ . » » • . \/ 
■ A ^\a) ii^ 2i^2\c j '^ '^ 

The h^s and ^'s are unknown coefficients, to be determined. 
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We have now to find the potential at any point in space. 
The mass of the spheroid o is jira^, and its potential is ^ira^/r. 

The potential due to the departure from sphericity, represented by the term in hi in 
the first of (17), is 

"~W" 2?^ \c / \r/ r*' ....... (18) 

This is written in a form convenient for passing to the case of r = a. It may also 
be written in the form 

"^ \c/ 2f^2W^-^' ....... (19) 

when it is in a suitable form for application of the transference formula (4). 
We shall now introduce two new symbols, namely, 

_ faV _ /AY 

Then (19) may be written 

and, of course, the similar potential with the other origin is 

3^^ [cj 2i^2 M^^^^' ....... (21) 

The whole potential at any point of space consists of the potentials of the two 
spheres and of the inequalities on each. The potential of the inequalities of the sphere 
may be written in the form (18), and of sphere in the form (21). 

Thus the whole potential is 



^-«^-^+-3r.!.2^2 c t 'r^ ..... (22-.) 



AnrA^ e ^Tra^ /A^ 'ZT SRiT^-^ &Wi ^ / . - \ 

'^ dc ' B'^ Z \cj ,?2 2^-2 li^^^^^ ..... (22-u.) 

The first line of (22) refers to origin o, the second to origin O, and to this latter 
half the transference formula (4) must be applied. 

Now apply (2) to the first term of the second Hne, and (4) to one term of the series 
in the second term, and we have 

AnrA^ G 4iitA^ ^'r,'^ /^Y Wk 

AC R 3c !c = o \oJ a* 

and 



cj 2^ ~ 2 



^2i + i— 3g (J 2^-2 tL h\i\ \c) a" 
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Thus the second line of (22) when transferred is 



47r-4 
3c 



■* = °° /aV m , s_ MS 'Zr ^^'" k + i\ P-^ faV w^' 



( Zi4~~'llt ) 



Then (22~i.) and (22~ii.) together constitute the potential now entirely referred to 
origin o. 

We want to choose the ^'s and JJ's so that each spheroid may be a level surface, 
save as to the outstanding term of the first degree. 

In order that (17) may be a level surface, when we substitute for r its value (17) in 
(22), the whole potential must be constant. In effecting this substitution, we may put 
r = a in the small terms, but in f ira^/r we must give it the full value (17). 

The constancy of the potential is secured by making the coefficient of each harmonic 
term vanish separately — excepting the first harmonic, which remains outstanding by 
supposition. 

We may consider each harmonic term by itself 

As far as concerns the term involving w^, we have, from (22-i.) and (22~ii.), as the 
value of the potential, 



a . 4:7rA 



3 



4 y-j^ I 

^ r oc 






and the value of r which must make this constant is 

but in the small terms inside [ ] we may put r = a simply. 

Make, therefore, the substitution, and equate the coefficient of w^ to zero. On 

dividing that coefficient by —~ — • f - j , we find 



' h T J- -— -L. 1 -X. .2- f _ \ V 



. _ ^k -\- "oT n +1 + 2(7) S 7. 1 -• 1 :: 7 H^i = 0- 

Inereiore 



h=i+^i^-j t^-jYiTi'^i-^'' ...... (23) 

and, by symmetry, 

\c/ i = 2 r.U! 4 -1 ^ ' 

Multiplying both sides of (24) by the coefficient of iZ", in (23), we have 



2Vc/ k\r\ r-1 ''~ ^[cl klr\r-l 



. (3^2 /«Y I^Y T " + ^"'+^' -^i^;.. 
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r = CO 



Performing S on both sides, and substituting from (23), 



r = 2 



h 



1 = i (^yrti^'^ + (If i^m: t^^m^t^kB^'- (-> 



0j\^2. ^1^1 ^ 



Now let 



(^,r) 

[^, % r] 



^ = 2 ^Irlr — 1 

% — L_ 

y = 2 T\i\h\r\ T -—X 



And (25) may be written 



J 



7i,= l + 



..>,r) + {f)^(-:Y(fY;j;[^,ir]^^/.. 



(26) 



• • (27) 



By imparting to A: all integral values from 2 upwards, we get a system of linear 
equations for the determination of the lus, and it will appear below that as many of 
them may be found numerically as may be desired. 

We now have to consider the series (26). 

Let 



and denote the operations 



^=1 



7 



B = 



i-r 



1 d)^ 



1 dh 



Wdy^ 



or 



XldV^ 



r^ by F^, 



Consider the function y^^E^, y log (1 + /3). 

Now 



Therefore 



Thus 



log (1 + /3) 



7' = 00 ^^^ - 1 

log(l-y)= S ^ 



r = 2 



r — 1 



1 ilk T— CO r<jk-\-r 



r = 00 



h ! dfy^y = 2 T —-1 



= t 



k -\~ rl Y 



r = 2 



kl rl r — 1 



(^,y) = i^^ylog(l + /8), (^,r) 



.s*.riog(i + £) 



« • 



(28) 



Next consider the function y ^E^EKylog (1 + iS). 
As before. 



and 



WE' vloa: (1 + iS) = " ^ t ^^^^ ^ — 

J^J^.yw^K^TP) ],^dr/,% i\r\r-l 



•)■• rr 00 



= % 



,=2 il r\ rl hi r 
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Hence 

{h, i, y] = - WM . 7 log ( 1 + ^), \h, i, r] = ^ B'E' . r log (1 + B). (29) 

We must now develop the symbolical sums of the series in (28) and (29). 
The following theorems are obvious :~ 

d"" f^ X p + ^ — 1 1 , , 

Then, by their aid, we have from Leibnitz's theorem — 

- /-Uog (1 + ^) = S^ ^-^-^-^ ~ y^- ^-^ log (1 + ^), 






— fsA^-i 



in which we interpret (— 1) !/(l — y)^ as log (1 + ^). 
Thus 

^^' '^^'^BoQc-t + i){k-t)t\k- ti ^'~'' ^^^) 

with/8°/0 = log(l+,8). 
Again 



Hence 



iltZotlk-tlk'-'t + 11 df{l ^ yy-^ 
i^o ^^0 1 1 k - tl rli -- rl k - t + 11 k -- t --' Hi + k - r - t + 



Ik h y] 

t=^'r=i 1 k + Hi + k-t+ 11 gi^Jc-r-t 

n: X X ■ — — — ' — — • . / ,^ 1 \ 

^^Q^^Q{iJ^k-T-t■i-l){i + k'-r-t) k-^t + ll t\k -t\ r\ i -- t\ ^ ^ 

In (30) and (31) the infinite series are replaced by finite series. 

From the form of the series it is obvious that the result must be symmetrica] with 
respect to k and % so that [^, % y] = [i, k, y]> ^^^ ^^^^ ^^ ^^^ obvious on the face of 
the formula (31). 

We shall find, therefore, the symmetrical form of (31) for the first few terms. 

If ^ = ^, r = i, we obviously have 

First term = {k + 1) {i + l)log (i + fi). 

The second term arises from t = k, r = ^— 1, and t = k-—l, r = i. The two 
corresponding values of (31) will be found to add together, and we get 

Second term = ^(k+ 1) (^ + 1) [2 {i + k) + Ik^ jS. 

3 D 2 
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The tliird term arises from t:=h, r-=i — 2; t =■ Tc — 1, r ■= i — 1 ', t=-h — 2, 
r = i, and we find— 

Third term = ~{h-\-l) {i + 1) 1^-^^^^ + * (^ - 1) 



3.2 



2! 



ih (* + 2) (A + 2) ik{i+k + l). „ 



2! 3! 



212 



_ A symmetrical form for further terms may be obtained by writing (31) first with 
i before k and then with h before % and taking half the sum of the two results. In 
computing these coefficients it is a useful check to compute from both unsymmetrical 
forms, when the identity of results verifies the computation. 

The following Tables have been computed from (30) and (31). The numbers are the 
coefficients of the quantities at the heads of the columns for the values of h and i 
written in the first column. The series (k, y) is terminable with jS^, and the series 
[k, % y] is terminable with jS^'^K 

In [k, ^, y] the coefficients have only been computed as far as ^^, so that the last 
which is given completely is [2, 4, y] ; however, with such values of ^ as we require, 
the series are carried far enough to give numerical results with sufficient accuracy. 







Table of {k, 


r)- 






Jc = 2 


Log (1 + /3) 


+ /3 


+ iS2 


+ /33 


+ /3^ 


+ /3'i 


3 


3 


1 

2 


• • 


• • 


* * 


h = S 


4 


6 


2 


1 

3 


• • 


• * 


h = 4i 


5 


10 


h 


1^ 


1 

4 


• • 


-k-h 


6 


15 


10 


5 


i| 


1 

5 



Table of \k^ % 7]. 



h=:2,i = 2 


Log(l + ^) 


+ /3 


+ ^2 


+ /38 


+ /34 


+ p^ 


+ ^« 


+ (37 


9 


27 


18i 


8 


14 


• « 


• • 




h=z2,i=^3 


12 


48 


46 


31 


12 


2 


• * 


• * 


Jc-2,i-4i 


15 


75 


92i 


85 


50i 


17 


2i 


» • 


7c=z%i = 6 


18 


108 


163 


190 


151| 


77f 


23 


&G. 


h^B,i = 3 


16 


84 


108 


103 


63 


22 


3i 


■ • 


^ = 3, f = 4 


20 


130 


210 


260 


219 


118 


361 


&C. 


h=zS,i=:6 


24 


186 


362 


552 


594 


434| 


206 


&C. 


h=z4i, i = 4i 


25 


200 


400 


625 


687i 


514 


248i 


&c. 


^ = 4, ^ — 5 


80 


285 


680 


1285 


1750|. 


1681 


1110 


&c. 


7<; — 5, * — 

, ,. 


36 


405 


1145 


2585 


4272 


50981- 


4345 


&'0. 
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We must now go back and determine the value of the outstanding potential of the 
first degree of harmonics, which will be annulled when rotation is imposed on the 
system. The potential is given in (22~i.) and (22-ii.); (22-i.) contributes nothing, and 
(22~ii.) gives us, for ^ = 1, 



3c 






w. 



a 



Thus, if we call Uy Ui the outstanding potential of the first degree, when referred to 
the two origins respectively, we have 



U^ 



U, 



4:7rA^ 
3c 

■~37 



i + %{^'^i±lr-^H?'-''^ -^ 



i = 9. i — 1 



i-2 



^ '' J ^a^'-l^ ' 



G 

c 



a 
W, 



y 



(32) 



5 3. The Potential due to Rotation, 



Intermediate between the two origins o and take a third Q^ and take the axes of 
l^and y) parallel to those of x and y^ and that of t, identical with that of z. Let Qo = d, 
QO = D. 

Then suppose that the system of the two spheroids is in uniform rotation about the 
axes of ^ with an angular velocity co. 

The potential fl of the centrifugal forces is given by 



But 



Hence 



fl = Ico^ {rf + 0. 









•^ 



X 



-3r = ■— fe 



• • • • • iooj 



. . . . . (34) 



XI 



1 ^3 ^2^3 ^« ^2 _ 2zd + d^) 



2"^ 2 



y^) + f (^^ + ^^ + ^^) "~ ^^d + d^\ 



Then, remembering that 



w^-= %' 



and if we put 
we have 



% 



X 



3 



1 /y»" ___ 1 nfii 

2*^ 2^ > 



^^ 



te^i 



^. 






f 2 = — ^0)^22 + 6^^% "~ co^d^te;^ + |a)V^ + \(i^d^. 



(35) 



Similarly the rotation potential, when developed with reference to the other 

origin O, is 

O = - ^0)2^3 + io)^ TTg - co^D TTi + ico^i?^ + Ico^Dl . . . (36) 



390 



PEOIESSOR G. H. DARWIN OK FIGURES OP 



The last terms of (35) and (36) are constants, and the term in r^, and that in W 
are symmetrical about each origin, and so the corresponding forces can produce no 
departure from sphericity in either mass; thus these terms maybe dropped. Next 
we have in (35) and (36) the outstanding potentials — o)^dw^ and — co^DW^, which 
will be annulled by other similar terms, and so need not be considered now. We are 
left, therefore, with the terms in q^ and Wc^, or in Q.^ and TF^. The g'^ ^^ ^ sectorial 
harmonic, the Wc^ a zonal, and it will be convenient to treat them separately. We 
shall begin with the zonal term. 



I 4. Disturbance due to the Zonal Harmonic Rotational Term. 

The potential whose effects we are to consider is ^cty^Wc^ or \oy^W<^, according to the 
origin which we are considering. 

If an isolated spheroid of fluid of unit density be rotating with angular velocity w, 
the ellipticity of the spheroid is Iboi^jlQir ; therefore we put 



15(1)^ 



e = 



16 



TT 



(37) 



Let us assume, for the equations to the two spheroids, 



r 
a 



3^ ^3 



i = 2 2^ — 2 \c 



I . -I -o-fc ' 1 ' . 



Jui 






S 



. . . (38) 



where k, Li^ are unknown coeflScients which are to be determined. We now have to 
determine the potentials at any point of the inequalities (38) on the two spheroids. 
The potential of the inequality |-€ Wc^jr^ in the first of (38) is 



4 3 3 1 3 



Wa^l-f"^. 



(39) 



The similar inequality in the second of (38) gives us 



W, 47r^3 ij\^ ^nj^ 



gTTZl . ^ . 36 ^g — 



Be 



B^ 



(40) 



The term in 4 in the first of (38) gives us, as in § 2, 



'37 2k - 2\g 









• • • • • 



(41) 
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The term in Li in the second of (38) gives us, as in | 2, 

^^^ \gJ 2i — 2 J?^^+i ^ ^ 

The potential due to rotation is ^co^t^^ or ^ox^W^, being the second term of (35) or 
(36) ; this term we find it convenient to write 

i^^^^Q'j- (43) 

The sums of the several terms (39), (40), (41), (42), and (43) are to be regarded as 
the potential of perturbing forces by which the spheroid a or the spheroid A is dis- 
turbed, and the arbitrary constants I and L are to be so chosen that each may be a 
figure of equilibrium. 

We may consider the spheroid a by itself, and the solution for it will afford us the 
solution for the spheroid A by symmetry. In order to find the disturbance, the formulae 
(40) and (42) must be transferred. 

Now by (4), with %■= 2, 



do ^^\g) JR' "" 3g ^^[gJ ^lo 2lkl \cj w 



fMtI --^=-^7-tMT >- "^TTT ;h;i- • • • (400 



And again, by (4), 

3c \g1 2i^2 E^'+^'~^ 3g ^\g) ,^o1T^" ^"^^ \J ^^' * ^^^ ^ 

Then (39), (40'), the sum of (41) from h= co to ^ = 2, the sum of (42') from ^ = 00 
to i = 2, and (43) together constitute the disturbing potential, all now referred to the 
origin 0. 

In order to find the disturbance of the spheroid a, we add the perturbing potential 
to iTra^/r, give r its value (38) in this term, put r = a in the perturbing potential, and 
make the whole potential constant by equating to zero the coefficient of each 
harmonic term. 

We will begin by putting r/a = 1 + 3^ Wc^/r^, and considering only the perturbing 
potentials (39) and (43). We have then, for the coefficient of w^/r^, 

Now, with the value of e in (37), 

— i^ra^' ^€= -— -r^oy^o? and — y% + i + i = 0. 

Hence the coefficient of %/r^ vanishes, and the term in e in (38) has been properly 
chosen to satisfy the perturbing potentials (39) and (43). 
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Following the similar process with the remaining terms of (38), and equating to 
zero the coefficient of lOjc^ we have from (40')5 (41)5 and (42% 



2Z; + 1 Uj, 



u "~" i!i 



/itC — - Ji 



+i.i^y^.±^'+if^r*r^'.!^i,=o; 



h\2\ 



G 1^2 ^- ^- *' — 1 



whence 






• (4^) 



By symmetry, the condition that the spheroid A may be a level surface is 






s. 



r\ 2! 



c 



2 t\ i\ i — 1 



t'^ • • • • * 



(45) 



Multiply both sides of (45) by f (-) ^^ - — r , and perform % on the whole, and 
substitute from (44) ; and we have 



I 



h 



if 



'A^h-V 2! 



' 



0) h\%\ 



= 5^ 






2 ! ^' ! ?M A ! ?" 



+ 






. . . . (46) 



Introducing the notation (26) for the series involved in (46), we have 



I 



k 



1 



'^Y, 



x(,+ i)(, + ,) + .|(|ygCM,r: 



+ (f) 



3 



a\3 /^\3^=^ ^, . _^ ry^-l 



C 



;. = 2 x/ — 1 



• • 



i = 2 



. . (47) 



Each value of h gives a similar equation^ and there is a similar series of equations 
with small and large letters interchanged. 
Now put 

^^ = -iV(7)'(^+i)(^ + 2)x, 



a 



3 



-^.• = iV: (^ + i)(^ + 2)A 



h 



and (47) becomes 



X^ ^ — 1 -|~ 



^ 



,• • • • • » • l4rOI 



a\^ /a\^ 



(k + 1) Qc + 2) \A 



)i(;)'ft2,r] 



+ (f) 



3/ „ 



aY /AY%^ (i + 1) (i + 2) r-, , ^-^ 7^-1 



c/ i^2 (7^ + 1) (7^ + 2) 



[k, i, r] -; r X,-. 



i-1 



• • 



• (49) 
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We attribute to h in (49) all values from oo to 2, and thus find a series of equations 
for the Vs. A similar series of equations holds for the A's. 

We must now find the outstanding potential of the first degree of harmonics. No 
such term exists in (39), (41), (43), but it arises entirely out of (40^) and (42'). If we 
write Vi for the outstanding potential, we have clearly 



i\ = 






c) 2! 11^7 '^^[c] -i^ i\l\i-l ' c a 



5^ 



whence 



% = 



47r^3 
3c 



10 ^ 



'a\^ 



CO ^ = 2 I — i 



A. 



a w. 



c a 



. . (50) 



and, by symmetry, 

^1 Q. 10 



3c 



-i-V € I -• 



A\2 

G 



\^/ i=2 '^ I 



c A 



. . (51) 



5. Disturbance due to the Sectorial Harmonic Rotational Term, 



In (35) and (36) we have found this term to be — ^ox^qc^ or ■— ^oy^Q^* 

We have already observed that, if the operation "T~i — t^ or S^ be performed on Wi^ 



the result vanishes when ^ = 1, 2, 3. 
Now, by (6), 



%D 



4> 



= t 



4! 



4' „ 4' 



1!^2! 



2!20! 



Hence \d?wjdp^ = 3, and, since S^^^,^ = i (a;3 — 3/2) cP'wJdp^, it follows that 



^3 = a;^ - 2/^ = |SH and ^, = 18^ W, 



(52) 



Hence the sectorial rotational term is •— •-^<^o)^hhv^ or — iVcy^S^TF^; this potential 
is of the second order of sectorial harmonics. 

Now, with € as defined in (37), let us assume as the equations to the two surfaces, 



T 

a 

E 
A 



= 1 



1 •- 



1 ^^ 



^ = 2 2^-2 



2'2; + iAA^* + i S2^. 



77?. 



^ + 2 



4 



6 



M^ 



A 



i = 2, 2^ — 2 



If. 






> . 



« + 3 



^' 



• (53) 



We have now to determine the potentials of the inequalities on the two spheroids 
expressed by (53). 
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The potential of the inequality -— y^-e Bhvjr^ in the first of (53) is 



K '-' /po " \ '}'' / '3'* 



---^•^^•i^-::j- = ~i^'^' :h:5^- (54) 



The potential of the similar inequality in the second of (53) is 



3g \c ""^"^"^ B^ ~ ' 



(55) 



The term in m/^ in the first of (53) gives us 

4cwA^ Sm^ faV faV "^ ^ S^'^^ 



+ 2 



3c 2^ - 2\c/ \r r^ 



. . . . . (56) 



The term in Mi in the second of (53) gives us 






Lastly, the sectorial term itself is 



12-0) a (■*-) ~^ ' ......,♦ (58) 



CCf T 



The sums of the several terms (54), (55), (56), (57), and (58) are to be regarded as 
the potential of perturbing forces by which the spheroid a, or the spheroid A^ is 
disturbed, and the arbitrary constants m, M^ are to be so chosen that they may each 
be figures of equilibrium. We may consider the spheroid a by itself, and the solution 
for it "will afford the solution for the spheroid A by symmetry. In order to find the 
disturbance, the formulae (55) and (57) must be transferred. For this purpose we 
require the second transference formulae. 

By (16), with % = 2, we have for (55) 



And by (16) we have for (57) 



3c I J 2^-2 i^2^+i 



3c M J ;, = 2 ^* - 2 ! /i; + 2 U - 1 '' [cj # 



Then (54), (55'), the sum of (56) from k= oo to k=2, the sum of (57') from 
i =z Qo to i = 2, and (58) together constitute the disturbing potential, all now referred 
to the origin o. 
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In order to find the disturbance of the spheroid a, we add the perturbing potential 
to f TraYr, give r its value (53) in this term, put r = a in the perturbing potential, 
and make the whole potential constant by equating to zero the coefficients of each 
harmonic term. 

*We will begin by putting r/a = 1 — -^-e h^wjr^, and considering only the perturbing 
potentials (54) and (58). We have then, for the coefficient of S^w^^/r^, 

Now, with the value of c in (37), 

iwa^ , ie — ^^o>%% and /4-- 1-^2 = 0. 

Hence the coefficient of B^wjr^ vanishes, and the term e in (53) has been properly 
chosen to satisfy the perturbing potentials (54) and (58). Following the similar 
process with the remaining terms of (53), and equating to zero the coefficient of 
§^'^^^+2? "^^ have, from (55'), (56), (57'), 



2^ + 1 

2^-2 



m 



Sm^ 



Ar6 



10^ 



c 



3 /?\^%°° 
2" 1 ) ^ 



h -^ i\ 



pi_i 



^^2^-2!/^^' + 2! i-1 



if; = o. 



or 



mj, 



'A\^ 



10 






^\ 3 i = 00 
2" ( "^ 1 ^ T 



h -\- i\ 



pi-i 



i = 2^ — 2!Z; + 2! i — 1 



■IXJ-i * 



(59) 



By symmetry the condition that the spheroid A may be a level surface is 



M, 



r 



10' 



a 



e 



+ t " S TT-. — r-. -: T mi (60) 



Multiply both sides of (60) by f (-) 



k + rl 



■pr-l 



r = 00 



^hJr^2lkT2l r -1' ^^^ perform^ S on the 



whole, and substitute from (59), and we have 



^k 



10^ 



'AV 



a 



2 



a 



3 >'= 00 



h -{- t\ 



Yr-l 



"^^^[c] Mc/ r^2 r-2\k + 2\ r-1 



T '\' i\h -^ r\ 



"*" ^^^ vJ \ J .fg ^^i ^2\r + 2\T -2\h + 2\ r 



pr~l 



7^ 



1 i — 1 



mu 



Now let us write 



y = 00 



{h, T}=i t 



h + t\ 



■\t — 1 



^2^ -Ilk -^21 



Kf^ 1j L 



r + il k -{• rl 



r = 00 

t 

r =.2 i - II r + 21 r - 11 k + 21 



V— 1 



SO that (61) may be written 



(61) 



(62) 



'A\^ 
c 



a 



% 



a 



3 



w/.= iV(t) +iV(5) l(r) {h v} + m 



'a\^ fA\^^ 



00 



t IM^I /"■%.'. (63) 



i = 2 



3 B 2 



396 



PEOFESSOR Q. H. DARWIlSr ON FIGURES OF 



Next put 



and (63) becomes 



'AY 



Jf; 



k 



10^ 



a 



s 



M,, 



e • 



. . . (64) 



i''=^+nf){^)\hT}+(w(^j{j)y\ij^Y-^,,. . . (65) 



We attribute to k in (65) all values from oo to 2, and thus find a series of equations 
for the jjLS. A similar series of equations holds for the M's. 
We now have to sum the series (62). 
Consider the function 

1 [(1 + ^)1^.^1]. 



k + 2 

' [(i+i8)^^^-i]=rf5[(i-r)-^~^-i]= ^ 



h + 2 



^4-2 



1 y = 00 
t 



Ic + 2 



k+2 , k+ 2 .k + S 



Yry + 



2! 



y + • • • 



k + rl 



Hence 



k+ 2r = 2 'k + llr — 1! 






{^.r}=T4^[(i+^r^-i] 



Next 



1 d}-^ 



{k+ 2)f i-l!^7^~3 



k+2 



wi{i+pr'-i\] 



(66) 



1 c?*-2 ^ = ^ 7^ + r! 



/+f 



7^^ — 1 1 ^7^-2 ,, = 2 ^ + 2!/' — 1 



1 r " 00 



k + r\ i + t\ 



-f r = 2. ^ - 11* + 2!r- llrH- 2 



17 



r + 3 



Henc( 



i-, i 7 



1 di-^ 



{k + 2)rf i-l\dr^i-^ 



{y^^[(i+i8r^-i]}. 



• (67) 



The differential in (67) must now be evaluated. We have, by Leibnitz's theorem, 



drf-^ 



{y«[(l+/3)^+^-l]} 



2 + 1 4^ 5^ 



^•-2! c?Y"^^ ^^""""^ 



^ + 1! 



J' = i — 2 






_ r! ^ — ?' — 2 ! dy ^7 



i-r-2 \ y) 



■y&-3 



v^ 4- S - 



i + 1\ i + k '— r — 11 7 



,i — f 4- 1 



&+1! (I-t)'**-*- 



■i + 1! 



-r + 



7'^ 



r = i - 2 

(1 _ yy~^ „? • • • '^ 



i — r--9> 



Substituting in (67), we have 



h i 7 



6^TT) 



r — i — 2 



i —- 2li + k -- Q" —- 1\ 



_i+(i + ^)... s„3! ^,.:^:-^, ::^^^ ^,--^,, ^ 



'« — r ■— 3 



r = 



.(68) 
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The following Tables, computed from (66) and (68), give the values of {k, y] and 
k, i, y I as far as ^ = 5, and k =■ 5, i = 5. 



Table of {k,y}. 



{2,7} 
{3,7} 
{4,7} 
{5,7} 



|[(l+/3)* 
|[(l+/3)^ 



-1]. 



Tahle of 1 k, i, y 



h 

h 
h 
k 

h 
h 
h 
h 



2, i 
2,i 
2,i 
2,i 

d,i 
3,i 
S,i 



2;i[(l+/3)* 
3;i[(l + yS)5 

4;t[(H-/3)^ 
5;f[(l + /3)7 






2;i[(H-/3)^-l]. 

3; f [(1 + ^)^1+ 1/3) -1]. 

4;f[(i + ;S)nn-|/3)-i]. 

5; [(l + ^)7(l + f/3)-l]. 



k 

k 

k: 

k 

k 
k 
k 
k: 



4, i 
4,i: 
4:,i 

5« 

:5, ^ 

5, ^ 

5* 
^ '2' 



2; i [(1 + ^)^-1]. 

3; i[(l + /3)^ (1 + 1/3)-!]. 

4; fra + /3)'(l + 3^ + fi^^) 

5; |[(i + mi + l/3 + -W') 



1]. 
1]. 



2; l-[(l+^)7--l]. 

3; f[(l + ^)7(l + f^)-l3. 

4;if[(i + mi + f^ + ¥/3^)-i]- 

5; f[(l + /3y(l + -V-^ + -V-/3^ + -V-^^) 



-1]. 



§ 6. Determination of the Angular Velocity of the System. 

The angular velocity of the system must now be determined in such a way as to 
annul the outstanding potential of the first degree of harmonics. 

Referring to origin o, we have from (35) — cw^ dwi directly from the rotation 
potential ; the remaining terms are Ui + Vi, since the sectorial harmonic term does 
not contribute anything. 

Thus, taking le^ from (32), and Vi from (50), we get for the potential 



0) 



3 



dwi + 



3c 7 






i = 2 



1 M^ 



+ io"^l ^ 



/ i=:2 '2^ — 1 \ 



>• 



Equating this to zero, 



3gA^ 
47r 



cl = ^M 



-*• "1 Q I 1 -^ . - A JLLi 



"T" 10 ^i 



MS 



l^ 






3 






>* 



(69) 
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And, by symmetry, 



i-rr 



D 



a 



N 



1 + 






lY'^h 



+^'(ij 



6 



a 



2 



H-t(fT|:^^^B^V-x.]}- (ro) 



Add (70) to (69), note that d + D = c, and -^^e = 3o)^/827r, and solve for o)^, and 
we have 



30)^ 

47r 



''Hf 



1 + 



^3(^8 



^ (A^ + a^) (^ i^2 i - 1 



'ti^iV^-'Hi + ry^-^Jid 



-Jl-1 



■f+C" 



3 
16 



'^^VaV -£(i±i)!(i±!)(,p-iA. + r,-^x.) 



(71) 



i = 2 



^ 



Now let 1 + -^ denote the factor by which {A/cf + i^l^Y i^ multiplied in (71). 
Then, if the two masses were particles, K would be zero, and (71) would simply be 
the usual formula connecting masses, mean motion, and mean distance in a circular 
orbit. Hence 1 -{' K is an augmenting factor by which the value of the square of 
the angular velocity must be multiplied if it be derived from the law of the periodic 
time of two particles revolving about one another. K, in fact, gives the correction to 
Kepler's law for the non-sphericity of the masses. 

This completes the solution of the problem, for we have determined the angular 
velocity in such a way as to justify the neglect of the harmonic terms of the first 
degree in §§ 2 and 4. 



§ 7. Solution of the Problem. 

We may now collect from the preceding paragraphs the complete solution of the 
problem. 

In (38) and (53) we have found that there are terms in r/a as follows : — 



1.^3 



1 S%^4 



3^ ^3 6*= ^3 



Now 



hence 



w^ = z^ 



L/y*'/C/ ^^^^ y, ni^ 



2 



■X 



y^, and h^w^ = 3 (cc^ — y^) ; 



2J/ ? 



^^^ _ I 8^^^ -%^- %x^ + ^^ = r^ - 3^^ 



and these terms are therefore equal to € (^ — x^jr"^). 

We note that € = 15a)7l67r = |a>V|7r, and that ofij^it is the ratio generally 
written m in works on the figure of the Earth. Then, from (17), (38), (53), the 
equations to the two surfaces are 
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X- 



3> 



a 



1 + € M — 7; •+ 



< j_\ 3 i = 002^4.1/^^ + 1 



r 



2 



Wi 



l + 4i 






^\3i=^co2^-+j. MV+1 

2U 



(^H-4);r-m 



S^'ii?,' 



« + 3 



rt^t 






W' ^3 117'. 



i?^ 



i2^ 



^- (72) 



From (27), (49), (65), we see that /%, ^3 . . . ^^ . . ., X^, X3 . . . X^ . . ., /x^, /^3 . . . /x^ . . ., 
are to be found by solving the equations resulting from all values of h from 2 to 
infinity in the following : — 



^z. — 1 =: 



-^^^*,r) + (|).0Y(frytM-,r]a., 



2 



X^ — 1 — 



3 



^\^ a 



(^ + l)(^ + 2)\c / \^ 



ft 2, r] 



(. H- 1) (^ f 2) (^)^ (0 ) t) .?3 ^^' ^' ^^ "• 



\^l 






a\^ a 



2 



'c&\3 IA\^^''zJ^ 



l(;)i5)P.r} + «)^(^)(-fj,S *..-,r 



i-1 



7' V' 



^x, 



. (73) 



and symmetrical systems of equations for obtaining the H'b, A's, and M's. 

With the values found by the solution of these equations we then evaluate K by 
formula (71) ; and we have 



3®^ 

47r 



'A\^ la 
J + [c 



r](i + in (74) 



We are now enabled to find the Vs and m's by the formulae (48) and (64), viz., 



fA\^ 
h = To € (A + 1) (^ + 2) { jj X^ 



^ , • • • • • • \i9jf 



J 



and the symmetrical forms give us the Ts and if 's. 

Having thus evaluated all the auxiliary constants, (72) gives the solution of the 

problem. 

It is well known that | X Zoi^j^^n is the ellipticity of a single homogeneous mass of 
fluid rotating with angu.lar velocity cw. Hence the first terms of (72) simply denote 
the ellipticity due to rotation in each of the masses, as if the other did not exist. 
Now the rigorous solution for the form of equilibrium of a rotating mass of fluid is an 
ellipsoid of revolution with eccentricity sing, the value of g being given by the 
solution of 

£ = cot^ g [(3 + tan^ g) g - 3 tan g].^^ (76) 



^ See, for example, Thomsoh and Tait's * Hatmral PMlosopliy ' (3), § 771, witli / = tang. 
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Hence it will undoubtedly be more correct to construct the surface, of which the 
equation is (72), by regarding the part of r under the symbol 2 as the correction to 
the radius-vector of an ellipsoid of revolution with eccentricity determined by (T6), 
where (jt>^/2TT is found from (74). 



§ 8. Examples of the Solution. 

The principal object of the preceding investigation is to trace the forms of the two 
masses when they approach to close proximity ; we shall thus be able to determine 
the forms when they are on the point of coalescing into a single mass, and shall 
finally obtain at least an approximate figure of the single mass. For this purpose we 
require to push the approximation by spherical harmonic analysis as far as it will bear. 
We shall below endeavour to estimate the degree of departure from correctness 
involved by the use of this analysis. The results will, therefore, be worked out 
numerically for such values of c/a as bring the two masses close together, and it will 
appear that the largest value of c/a assumed for numerical solution is such that the 
surfaces cross; in this case the reality will be a single mass of a shape which it will 
be possible to draw with tolerable accuracy. 

The computations are facilitated if, instead of assuming c to be an exact multiple of 
a, we take c^ a multiple of a^ ; that is to say, we shall take l/y as an integer, and 
therefore 1//8 also an integer. 

"We shall in the first instance suppose the two masses to be equal. In the following 
examples, then, we have A •= a, T = y, B = ^, and the two masses assume the same 
shape. 

The computations will be carried through in detail in two cases, viz., when y8 = y, 
and when /3 = ^. The results will also be given for ^ = ^* 

When /3 = I, y = i ^M = 2-8284, and when /3 = i, y = h ^A^ = 2'449. Thus 
the distances of the centres apart are 2f and 2f of the mean radius respectively. 
The numerical details of the two computations may be stated pari pass4, and the 
numbers applying to /3 = } will be distinguished by being printed in small type. 

In the case of /3 = ■^, we have y = j, c/a = 2*6458 ; but only the final result is 
given, without intermediate details. 

The first step is to compute the values of the several series by means of the Tables 
in S5 2 and 5. 

The numerical results are as follows. 
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Table of [Ic, y). 



Table of {h, y}. 





^ = f 


^ = 1 


y^ = 2 


•839 


1-167 


A; = 3 


1-433 


2-012 


7^ = 4 


2-204 


3-125 


7^ = 5 


3-163 


4-536 





^ = f 


^ = * 


h-2 


•177 


•268 


k-S 


•190 


-298 


7^-4 


•205 


•331 


k = b 


•221 


•369 



Table of [k, i, y 



Table of k, i, y. 





^ = f 


P = i 


k~-2, i- 2 


5-460 


7-847 


h = 2, i=B 


9-494 


13-895 


k = 2, i = 4 


14-875 


22201 


k=2, i=b 


21-780 


33-190 


k - 3, i = 2 


9-494 


13-895 


A; - 3, ^ - 3 


16-667 


24-969 


i^ = 3, ^ 4 


26-384 


40-517 


/^ = 3, ^ = 6 


39-047 


61-574 


/i;-4, ^ = 2 


14-875 


22-201 


7^-4, i - 3 


26-384 


40-517 


7^-4, ^' - 4 


42-214 


66-840 


k — 4^, i — b 


63-183 


103-372 


k- 5, i- 2 


21-780 


33-190 


k~-5, i-3 


39-047 


61-574 


k — 5, i — 4i 


63-183 


103-372 


7i; rr 5, ^ = 5 


95-690 


162-831 





|3 = f 


^ = h 


k ==z 2, 


i=:2'\ 


-177 \ 


-268 


k == 2, 


^ = 3- 


-475 


:744.. , 


k =- 2, 


i = 4 • 


1-024 


1-655 


k == 2, 


t = 5 ; 


1-933 


3-229 


^ == 3, 


i = 2 


-190 


•298 ' 


7' ==3, 


' ^ = 3 


-519 


-844 ' 


k - 3, 


^ = 4 


1-137 


1-921 


^ == 3, 


i:= 5 


2-183 


3-837 


k-4, 


i-2 ; 


•205 


•331 


k--A, 


i: = 3 - 


-569 


-961 


k -'- 4, 


i = 4 


1-266 


2-238 


k=:^4<, 


^•=5 


2-471 


4-578. - 


k =: 5, 


^-2 


-221 


-369 . 


k -". 5, 


^• = 3 


-624 


1-096 


7i; =: 5, 


i 4 


1-412 


2-616 


^ = 5, 


^- 5 


2-803 


5-479 



Witli these values for the series, we have to compute the coefficients of the systems 
of simultaneous equations (73). The equations lend themselves more readily to 
solution if we consider ^^ — 1, X^- — 1, /i,- — 1, as the unknowns instead of hi, X,-. /x,-. 
The results are given in the following equations. 
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The upper coefficients correspond to the case of ^ 
small type, to the case of ^S = -|-. 



1 

7 



; the lower ones, printed in 



k 



2 



h 



li. 



1 = 






•05902 + -00800 (^^-l) + '00033 (^3-!) + -00004 (/^.i 

•13516 -01362 -00201 '00036 



K-1 



•10086 + -00522 

-23385 -02412 

•15529 + -00817 

-36167 -03854 

-22321 4- -01196 

•53150 -05762 



+ -00057 

-00361 

-f -00091 

•00586 

+ -00134 

•00891 



+ -00008 

•00065 

+ -00012 

•00108 

+ -00018 

•00165 



1) + -00001 (li^ - 1) + 

•00007 



+ -00001 

•00012 

-f '00002 

•00021 

+ -00003 

•00033 



+ 



+ 



+ 



X. -1 = 



Aq 



"4 



1 = 



1 = 



X.- 1 = 



•06400 + -00300 (X3-I) 4- -00054 (X3-I) 4- '00011 (X^ 

•15158 ^01362 -00335 -00089 

•06677 + -00313 

•16114 -01447 

-06976 + -00327 

•17174 ^01542 

•07297 + -00342 

•18352 -01646 



+ 


•00057 




•00361 


+ 


•00060 




•00391 


-- 


•00064 




•00424 



+ 


•00011 




•00098 


+ 


•00012 




•00108 




-00013 




•00118 



1) + ^00002 (X5 - 1) -{- 

-00023 



+ -00002 

•00026 

-F -00002 

•00029 

+ -00003 

•00033 



+ 



+ 



+ 



;^2 - 1 = -01184 + '00010 (^^3 - 1) + -00003 (/X3 - 1) + -00001 (/^^ - 1) + -00000 (/^^ - 1) -h . . 

•02818 -00047 ^00022 '00008 *00003 



^^ _ 1 - -01274 + '00010 


+ -00004 


+ -00001 


+ -00000 


+ ... 


•03127 -00052 


•00024 


•00009 


•00003 




/^^ 1 - -01374 + '00011 


+ -00004 


+ '00001 


+ -00000 


+ ... 


•03478 -00057 


•00028 


-00011 


•00004 




^g ~ 1 = -01482 + -00012 


+ -00004 


+ -OOOOI 


+ -00000 


+ ... 


•03879 -00064 


•00032 


•00013 


•00004 


> « . • 



The solutions of these equations are obviously found by an easy approximation ; 

they are 

X^=: 1-0642 
1-1544 



U^ — 


1-0593 




1-1377 


H- 


1-1012 




1^2382 


K 


: 1-1559 




1;3719 


K 


: 1-2-241 




1^6424 



11^— 1-0118 

1^0282 



X 



= 1-0670 

1^1642 

\— 1-0700 
1-1750 

X, = 1-0732 

1-1870 



/Xg 



/^4 



/^5 



1-0127 

1-0313 

1-0137 

1^0348 

1-0148 

1-0388 



the small figures corresponding, as before; to the case of ^ ~ 



5- 



EQTJILIEEIUM OF ROTATING MASSES OF FLUID. 



403 



With these values of the h's and X's, I find 



ili.i^iT,.-- .Q^iQ. 2v^i^ti}ili:ii).i\.-i.QQ.o. bV^^^ 



2t-~-~Y^^hi=^ -8718; 
^""^ 1-3005 



i — 1 



/X.=z 1-3949; , 

2-8542 ^ ^■ 



the summations, of course, stopping with i = 5. 
Applying these in (71), we have, when 



y8 = f, 1 + K = 
whence 



1 + -02891 
1- -00880 



= 1*0380; or, when ^ =:= I, 1 + ^ 



5^ 



3o)3 



•08839 X 1-0380 = -09175. 
•13608 X 1-0877 -1481 



•00829, 
•01701 



1 + -0664 



1 - -0195 



; = 1-0877, 



Thus the angular velocity of the system has been found. 
Next we have 



~^^e[~] = '001434. 
^^/' -00309 



Introducing this into (48) and (64) with the previously found values of the X'sand /xs, 



/g z= '0183, m^ 
•0428 

4 == -0306, mg 
•0720 

Z^^ =: -0460, m^ 

•1089 



h 



= -0646, m 
•1541 



•00145, 

•0032 

•00145, 

•0032 

•00145, 

•0032 

•00145, 

'0032 



r h + k 



> ; and hence < 



/^3 + 63 



K + h 



h + ^6 



v~ 



1-0776, 
1-1806 

1 • lolo, 
1-3100 

X iy u 1 y, 

1-4808 

r2887. 
1-6965 



By taking the differences of A + ^? we may conclude that 

/^g + Zg=l'39, 
1-96 

and this sixth harmonic term will now be included. 

It appears from the values of the w/s that the harmonics of the type S%tV+2 ^^^® 
practically negligible, excepting the term 8%^, and that in that we may neglect the 
part depending on mg. 

Now, if r denotes the radius-vector due to the rotation, and Sr the increase of 
radiuB-vector due to the mutual influence of the two masses, v/e have 



a 



IDo 



IC 



^•i , .AAO^% 



V), 



'1191 ^^ '0309 -f + •0100—'+ '0035 "f + -0013-^+ , . . 



•2008 



'0637 



•0252 



•0108 
3 F 2 



■0048 



. . (77) 
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We next have to consider the values of r, the radiuB-vector of the eUipsoidj due to 
rotation. 

We might compute from the spherical harmonic formula 

The results so computed will be compared with the others computed as shown below. 
The following Table of the angular velocity and corresponding eccentricity e of the 
equilibrium ellipsoid of revolution is extracted from Thomson and Tait's ' Natural 
Philosophy/ §772 : — 





27r 


•3 


■0243 


•4 


•0436 


•5 


•0690 


•6 


•1007 


■7 


•1387 


•8 


•1816 



From this we find by interpolation that, when 3o)^/4it = '09175, 6 = '472 ; and, 
when 3a>V47r = -1481, e = -594. 

These, then, are the eccentricities of the ellipsoids whose radins-vector is r in the 
two cases yS = y, ^ = -|^. 

The equations to the generating ellipses are 

r 1 - -0806 . ^ . . r 1 - -1353 



a 



1 - -2228 cos^ ^ for ^ - ^ , and ,, - ^ _ .3535 ^^^2 ^^ ^ - i 



The following are the computed values of rja for each 15"^ of 6, the latitude, the 
small figures written below appertaining to the case of yS = \, 

0= 0° 15'^ 30^ 45"^ 60^ 75'' 90° 

^ = |: ~ = 1-0429, 1-0330, 1*0074, -9753, '9461, '9264, '9194 

^ = 1: 1-075, 1-056, 1-009, '953, '906, -875, -865 

Computing from the spherical harmonic formula, I find 

^= f: ^==1-0382, 1-0305, 1*0096, '9809, ;9522, *9312, ;9235 

.p = -|: 1-0616, r0490, 1*0154, -9692 -9230, -8892, '8768 



60° 


75° 


90° 


- -031, 


— -055, 


— -056, 


- -060, 


- -094, 


- -092, 
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The greatest discrepancy occurs when /3 = I and 6 = 90°, and the difference 
between the two results is fjj of either. It foUow^Sj therefore, that in drawing the 
figures it is not of much importance which results we take. But, as above remarked, 
the radius-vectors computed from the true ellipsoidal figure are the more correct. 

The formula (77) for Sr consists of a series of zonal harmonics. The pole of 
symmetry of these harmonics lies in the equator of the elHpsoid of revolution defined 
by r, and is that point of each mass which lies nearest to the other. Then, denoting 
by co-latitude estimated from this pole, I find that the numerical values of Sr for 
each 15° of ^ are as follows : — 

0= 0° 15° 30° 45° 

/3=\: -~ = + -1G5, +-141, + -084, + 'OlO, 

ft= li + -280, + -257, + -142, 4- '024, 

105° 120° 135° 150° 165° 180° 

— -037, — -004, + -032, + -065, + -088, + '096 

-•059, --002, +-055, + -106, +'143, 4- 'ISo^ 

These have to be combined with r, so as to give the radius-vectors of the mass of 
fluid along two sections, one perpendicular to the axis of rotation {which may be 
called the equatorial section), the other through the axis and the two centres (which 
may be called the section through the prime meridian). Taking the case of ^ = y, we 
add the successive values of Sr to the equatorial value of r, viz., 1*043, and thus find 
the equally-spaced radius-vectors along the equatorial section. Next we add the 
successive values of Sr to the corresponding values of r, and thus find the equally- 
spaced radius-vectors along the prime meridian. The results are as follows : — 

(9= 0° 15° 30° 45° G0° 75° 90° 

Equator, — ^' = 1*208, M84, 1126, 1*06% 1*012, '988, '987, 
Pr. Merid. = 1*208, 1*174, 1*091, '994, *915, '871, '863, 

105° 120° 135° 150° 165° 180° 

r006, 1*039, 1-075, 1*108, 1*131, 1*139 

•889, *942, 1*008, 1*072, 1*121, 1*139 

These results apply to the case of /2 = y ; those for ^ = |- are found in the same 
way, and are given in the figures referred to below. 

When /8 = y the distance between the centres is given by c/a = 2*828. I have 
also worked out the case of /3 = |, although none of the numerical details are given 

here. 

In figs. 2, 3, 4, and 5 (Plates 22 and 23) are exhibited the figures which result 
from some of these computations. 
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Figs. 2 and 3 refer to /3 = ^-, 4 and 5 to that of yS = -^, and the numerical values for 
yS = y, given above, make it easy to draw a figure for fi = y. 

Since in these cases the masses are equal, the two halves of the figure are the 
images of one another. The numerical value of each radius-vector is entered on the 
plates ; and other numerical data and explanations are given. 

Figs. 2 and 3 correspond to yS = -|-, and here the figures as computed cross one 
another. The reality must, therefore, be two bulbs joined by a stalk, like a dumb- 
bell. The dotted lines have been filled in conjecturally, and must show pretty closely 
what that single figure, formed by the coalescence of the two masses, must be. 

Figs. 4 and 5 show in a similar manner the case of ^ = ^, and here the two masses 
are separate, although nearly in contact. When /3 = y the shapes present similar 
characters, but are wider apart. 



§ 9. On the Use of Spherical Harmonic Analysis as a Method of Approximation, 

Spherical harmonic analysis gives less accuracy as the bodies considered depart 
more and more from spheres. How far, then, do our results present an approach to 
accuracy ? To answer this question, we have to find how nearly the potentials at the 
surfaces of these figures may be computed from the spherical harmonic formulae. 

It would be laborious to make an accurate computation of the potential, and it 
fortunately appears to be unnecessary to do so, since a sufficient answer may be 
obtained in another way. 

The potential of an ellipsoid of revolution may be computed either rigorously or by 
harmonic analysis. With a certain degree of eccentricity the approximate result will 
agree badly with the rigorous one. 

If the ellipsoid consists of a fiiuid of unit density, there is a certain angular velocity 
which makes it a level surface. If o) be that angular velocity, then we know that 
the spherical harmonic solution would give 1 — ■ Iboi^jlQir as the ratio of the minor to 
the major axis. If then c, a, are the rigorous values of the minor and major semi-axes, 
the harmonic approximation is good if cja does not differ much from 1 — 15^^1677. 

If we denote by 1 — /x the factor by which the approximate value of the ratio of 
the axes is to be multiplied in order to obtain the rigorous value, we have 

^ 1 - l5co3/167r ' 

and /x may be regarded as a measure of inaccuracy. 

A table of the values of co^j^rr^ corresponding to various eccentricities 
c =z- v^^C^ — i'^'I^Y)^ ^^ computed from the transcendental equation in Thomson and 
Tait's ' Natural Philosophy,' § 772. From these I compute as follows :— 
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e = 


-('-S) 


c 

a 


^ ~ "l67r 


Difference. 


1 




•1 


•9949 


•9949 


•0000 


Large 




•2 


•9798 


•9799 


•0001 


9799 




•3 


•9539 


•9544 


•0005 


1909 




•4 


•9165 


•9182 


•0017 


540 




•5 


•8660 


•8705 


•0045 


193 




•6 


•8000 


•8111 


•0111 


73 




•7 


•7141 


•7399 


•0258 


29 




•8 


•6000 


•6595 


•0595 


111 




•9 


•4359 


•5869 


•1510 


3-9 



The measures of inaccuracy corresponding to the values of e in the first column, or 
the values of cja in the second, are the reciprocals of the numbers in the last column* 
We thus see that there is still a considerable degree of approximation when e= '8, 
or when the ratio of the axes is 3 to 5, for the measure of inaccuracy is y^ ; but for 
6 = '9 the approximation is bad. 

Now the shapes of certain egg-like bodies have been computed by the spherical 
harmonic method, and it seems safe to assume that the approximation has given 
about the same degree of accuracy as would hold in the case of atl ellipsoid of 
revolution whose minor axis bears to its major axis the same ratio as the shorter axis 
of the Qg^ to the longer. 

Turning now to our computation, and considering only the more elongated or 
meridional sections, we see that, when yS = y, the longer axis is 1*3 55 + 1*2 30 = 2*585, 
and the shorter 2 (1 — '227) = 1*546 ; and the ratio 1*546 : 2*585 is *6, which corre- 
sponds to the measure of inaccuracy 1/11*1. It might, however, be more legitimate 
to adopt two different measures, and at the pointed end of the egg to take the ratio 
•773 : 1*355 = *575 which will correspond to a measure of inaccuracy about xb-; and at 
the blunt end to take the ratio *773 : 1'230 = *63, which would correspond to a 
measure of inaccuracy Y2 ^^ A- 

In the case of yS = -|- the two masses cross one another, and the result has been 
used to give an approximate picture of the dumb-bell figure of equilibrium. We now 
see that even in this case there is a sufficient degree of approximation to give a very 
good idea of the accurate result. 

In the case of the meridional section, where yS == |, we have for the ratio of axes at 



the pointed end of the Qg^ 

1 



1 - -1717 

1-2712 
- -1717 



■8283 



2-5-; ^i^d at the blunt end ■^.^^^r^ 



1-2712 

•8283 



- zz: '65, and measure of inaccuracy about 



1-1746 



*7l, and measure of inaccuracy 2V 



•863 



1 — '137 
In the case of yS = y the similar figures are, for the pointed end, "i^o aq~ ^ 1 !o a q ^ '^^^ 



1-208 1-208 



and measure of inaccuracy about -^^ ; and for the blunt end 
measure of inaccuracy perhaps about y^-. 



137 -863 



1-139 



1-139 



-= ^^iS^ and 
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It thus appears that as the bodies recede the accuracy increases with great rapidity, 
and in the two cases considered last it is hardly necessary, from a physical point of 
view, to consider greater accuracy than that attained. 

It must be remarked, however, that this way of estimating the degree of inaccuracy 
must necessarily give much too unfavourable a view. 

If we have a single mass of fluid departing considerably from the spherical form, it 
is clear that the potential computed on the hypothesis of a layer of surface density on 
the true sphere will come to depart largely from the potential at the surface of the 
fluid. If, however, we compute the potential of such a mass at points a little remote 
from the surface, the approximation will be much closer. Now, where there are two 
masses, as in our problem, the potential at the surface of either mass consists of two 
parts, one due to the mass itself, the other due to the other mass. As regards the 
first of these two parts, the above criterion is applicable, but as regards the second 
Dart it gives too unfavourable a view. 

Now in the case of the single mass the deformative forces due to centrifugal force 
are considerably vitiated by computation at the spherical surface instead of the true 
surface, whilst in the case of the two masses the tide- generating forces are computed 
with greater accuracy than is shown by the criterion, 

Under these circumstances it has appeared worth while to give another figure 
below, which, judged by the criterion, would be no approximation at all. 

The reasons for giving this figure will be stated when we come to it. 

I 10. To find the Moment of Momentum of the System. 

Rotating figures of equilibrium are classified according to the amount of moment of 
tnomentum with which they are endued. It is, therefore, interesting to determine the 
moment of momentum of the systems now under consideration. 

We must begin by finding the moments of inertia of the two masses. Let S/, St, 

denote the moments of inertia of the shells of zero mass lying on the mean spheres of 

radii A^ a. 

Then 

8i == \^ (jj^ 4" z^) (r — a) a^ c^ct, 

where dz!T = mn 6d6 d(f), and where the integral is taken throughout ^.ngular space. 

Now 

2/^ + ^^ = |«^ + i«^ (5 - i^-^^ 

and r — a is the sum of a series of harmonics. Then, in consequence of the properties 
of harmonic functions, we need only consider the harmonics of the second degree in 
r — a, and 



Si = ^a^ 



3^ 






+ 1 



rt-.2 I ' -4 \ 2 ^»iJ / 
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But 






nr 



4 



5^5 



2" ^3 



C^OT 



Itt . 3, 



36 — 



Odd" 

IGtt 



and the moment of inertia of the mean sphere is --r^rra^ ; hence, if we write 



om^ 






i^ = 






^\3 



1(^2 + y + 3m J 



[1 (^, + i,) + 3.¥J 



the moments of inertia, i and /, are given by 

I = -^^tra^ {I -^r F). 



We already have in {71) 



3«2 

Air 



'fi + (: 



(l+iT). 



Hence the sum of the rotational momenta of the two masses is 



(i -|- Z) CO = 



2. 
5 



'4:7r\^ 



d 



[a^l+f)+A^l+F)'] 



aV (A 
e +17 



I 



{l+Kf. 



The whole system revolves orbitally about the centre of inertia with an angular 
velocity w : hence the orbital momentum is 



But 



d^ 



a^ + A^ 



Hence the orbital momentum is 



and this is equal to 



A^ 4- a^ 



D = 




arc 


a;' 


+ ^3 


(t)d^, 










It will be convenient to refer the mass to the radius of a sphere of the same mass 
as the sum of the two. 

Let h be the radius of such a sphere ; then 

W zzz A^ '\' al 



Thus the whole moment of momentum is 



(i^)^ h' f 



a' 
h 






(1+/)+ T {l+F) 



!)■+ fir (f )^ (f ' 



(1 + Kf, 



We shall therefore compute the coefficient of (Itt)' &l 
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Computing from this formula, I find the following values of the moment of 
momentum in the case where the masses are equal, when 

Now I find by a numerical investigation"'' that, if we imagine a mass of fluid equal 
to %7rW rotating in the form of a Jacobian eUipsoid of three unequal axes, then, when 
the momentum is {%7r)^h^ X '392, the axes of the ellipsoid are 1*8986, 0-81136, 0'6496 ; 
and, when the momentum is (fir)' 6^ X '644, the axes are 3*1366, 0*5866, 0*5456. 

It seems probable, then, that the Jacobian ellipsoid of mass ^irW becomes unstable, 
at least as soon as when the moment of momentum is somewhere about {-%Trf ¥ X '5. 

It may be worth mentioning that the greatest moment of momentum for which the 
ellipsoid (of mass ^ttW) is stable, when it is a figure of revolution, is (Itt)' 6^ X '3038. 



I 11. 0)1 the Conditions under luhich the two Masses may he close to one another. 

If at any point on the surface of either mass the sum of the tide-generating and 
centrifugal forces is greater than gravity, it is obvious that equilibrium cannot subsist. 
It is also clear that, if this condition is to be found anywhere, it will be at that point 
of tbe smaller mass "which lies nearest to the larger mass. Hence, in order that the 
system may be a possible one, we must satisfy ourselves that at that point gravity of 
the body itself exceeds the sum of the tide-generating and centrifugal forces. 

To determine the Umitations of size and proximity of the smaller of the two 
masses to a high degree of approximation would be very laborious, and we shall, 
therefore, content ourselves with a rough investigation, to be explained below. 

We shall now find approximations for the shapes of the two masses and for their 

potentials. 

The radius-vector of either mass and the potential may be expanded in powers of 
aje and Aje, and a term involving c'' in the denominator will be referred to as being 

of the n*^ order. 

Now the term of the highest order which can be included without the introduction 
of great complication is the 7th, and we shall content ourselves with that term. 

The expressions for the various parts of the potential have been developed above, 
but it may be observed that the terms involving the first order of harmonics may be 

* 'Roy. Soc. Proc.,' vol. 41, 1887, p. 319. 
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omitted, since' they are subsequently annulled by a proper choice of the angular 
velocity. 

From (22-i.) we have 



„ I, ij 

3 r """ 3c 2 






,> -J-' . . t 



The last term in the development to the 7th order is that involving tVQ, Then it is 
clear that we require ^2 correct to the 4th order, h^ to the 3rd, and so on. But (25) 
shows us that the A's are equal to unity to the 4th order inclusive. Hence, in the 
above, all the h^s may be treated as unity. 

Again (22-ii.) when written in reference to the origin o affords other terms, in 
which all those included under S2 are of the 8th and higher orders, and negligible ; 
and the rest (with omission of the first harmonic term) gives 



4.7rA' 
"IT 



a\^Wc 



a\3 w, 



3 



c 



a' 



+ (;)^ + --- 



Thus this first part of the potential is, to the 7th order inclusive, 



3 r 3c: h^^Ac. 



a 



k + i 



r 



k 



2{k-l)\:r 



+ - 

a 






Next, from the expression for O in § 3, we have a term in the potential due to 
rotation + ^wV^. The remaining terms due to rotation will be taken up later. 
From (71) we see that, to the 7th order inclusive, 



3 



ft) 



3 



'A\^ 



a 



\3 



A-nr 



+ <; 



Hence m^ and e are of the Si-d order ; and from (48) and (64) it follows that the 
factors by which the fs and ms are derived from the X's and ju,'s are of the 5th orden 
And, since the X's and jit s only differ from unity in terms of the 5th order, it follows 
that the Z's and m's are of the 5th order. Then (41) and (56) show us that all the 
terms in I and m are negligible. 

The first set of terms due to rotation and to the corresponding deformation are 
given in (39) and (43), and together contribute 



10 9 



3 
2 



ri \a, 






The second set of terms due to rotation, and to the corresponding deformation, are 
given in (54) and (58), and together contribute 



^iii^a?" 



3 
2 


3 G 2 


1 s%,i 
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Hence, to the 7th order inclusive, we have 



F = 



4z7ra^ a . 4:7rA^ ^ = ^ faV 



3 T 









2(k-l)\r) ^ W J ^'^ 



W 



3. 

2 



a 



A^ 



3 






2' 



a, 



VJo 



W'^^ 






. . . . (* ^/ 



Now the expression (72) for the radius-vector of the mass a to the same order of 
approximation gives us 



a 



I i^^ / ^g -^ ig%A /:4\^'|'2^±i/^Y"'^!^^, 



16 



TT 



mf^ 



G i^ = 2 2/<? — ^\^. 



11^ 



and a similar expression for RjA, 

To determine the inward force at the pole of the mass a, where it is nearest to the 
mass A, we must evaluate — dV/dr, and in the first term substitute the above 
expression for r, and in the remaining terms put r = a; also at this pole iv^/r^ = 1, 
and 8hv^ = 0. 

Then, difterentiating (78), 
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Substituting this for the first term of (79), we have 
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Thus the criterion of the possibility of equiUbrium is that 
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should be positive. 

Bat the radius-vectors of the poles are 
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Therefore 
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Now the interval between the two masses is c — {r + E) ; hence, if the two masses 
are just in contact, 

c=a + A+^ l{7a + A)A^ + {a + 7A)a'] + £, A\t' {A + a) 
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In order, then, to test whether equilibrium is still possible when the two masses are 
just in contact, it is necessary to determine c from (81) ; and then, substituting in (80), 
find whether C is positive or not. 

The solution of an equation 
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can only be performed by trial and error. 

Now suppose that the solution is Cq + Sc, where Sc is small ; and that 
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Then it is obvious that 

Bo 6'-^ — Cq 



^0 T^ . 3 + "4 ^ • • • 
^0 *^{) 



and 



8e=('^ + ^^+ . . .)-, and 0=:^ C,+ 8C 



^0 ^0 / ^0 



It is not hard to find an approximate solution Cq by trial and error, and the correct 
results may then be found thus. 

Consider the case where the two bodies are equal to one another, and put a = ^ = 1. 
The equations then become 

^ 3c^ ^ 2c' ^ c'^ 4c6 ^ 5^7 

14^_7_^^^11 18 

3c3 "~ 2^ "~ 2^5 "" 2^6 ■" 2c7 * 

By trial and error we find c = 2*535, (7= + '557. 

From this we conclude that equilibrium still subsists when the two masses are in 
contact. 

When a = J = 1, c = 2*535, we have y = (a/cY = r^a^ ^^^^ ^ = y/(J —y) = "s'Vs"- 
Our figures showed that when ^ = ^ the two masses were nearly in contact, and when 
/3 = -5- they crossed. 

This result is, therefore, in accordance with the figures. 

Next pass to the case of an infinitesimal satellite, and suppose a infinitely small 
compared with A and c, and that A = I. The equations are 



c — 
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43 
12c3 



The solution of the first equation is c ::^ 1 *226, and this value of c makes (7= — '94. 
Hence we conclude that an infinitesimal fluid satellite cannot revolve with its surface 
in contact with its planet. 

C vanishes when c == (ff)^ =^ 1*89. Hence it appears that the nearest approach of 
the infinitesimal satellite to the planet is 1*89 mean radii of the planet. The nature 
of the approximation adopted is, however, such that in reality the satellite must 
lie further from the planet than this, perhaps at two radii distance,"^ The satellite and 
planet of which we here speak are, of course, supposed to revolve as parts of a rigid 

* [See Roche, * Montpellier, Acad. 8ci, Mem./ vol. 1, 1847''"1850, p. 243 (added Oct. 5, 1887).] 



EQUILIBRIUM OF ROTATING MASSES OF FLUID. 415 

body. Now, if for equal masses equilibrium still subsists when the two masses are in 
contact, whilst for infinitesimal mass of one equilibrium is impossible with the masses 
in contact; it follows that for some ratio of masses equilibrium can just subsist when 
they are in contact. 

The question, therefore, remains to determme this limiting ratio of masses. 

I find, then, that when a = 1, A = 3'4, we have 

= 4-4 + [2-25684] c~H[l '94945] c-H [2-07156] c"-s + [2-3761 9] c-^ + [2-80737] c--^, 
0=1 — [2-15205]c-^-[2-13850]c-*^[2-24765]c-^ — [2-33480]c-^-[2-40735]c~^ 

the numbers in [ ] being the logarithms of the coefiicients. 
The solution of this is c = 5*57, which makes (7= — '006. 
Again, when a =: 1, ^ = 3'3, we have 

c = 4-3 +[2-21556] c~^+ [1-91353] c-^+ [2-03266] c-^ + [2-32731] C-6+ [2-74467] c-^ 
(7=1 — [2-11347] c-3- [2-09961] c^^~ [2-20876] c-s- [2-29591] c-6 -[2-36846] c-^ 

the solution of which is c = 5-45, which makes (7 = + '010. Since (3-4)^ = 3 9 '3, and 
(3-3)^ = 35-9, it follows that the ratio of the masses in the first case is 1 : 39-3, and in 
the second I : 35-9. 

From this it appears that when the ratio of the masses is about 1 to 38 equilibrium 
is still just possible when the two masses touch. 

It must be borne in mind, however, that the nature of the approximations adopted 
in this investigation is such that the results in this limiting case are only given very 
roughly, and it is certain that actually the limiting size of the smaller of the two 
masses must be greater than as thus computed. 

We can only conclude that the limiting case occurs when the ratio of the masses is 
about 1 to 30, or the ratio of the radii about 1 to 3. 

There is one other case which it is interesting to consider, namely, to find the 
limiting proximity of the Moon to the Earth, both bodies being treated as homogeneous 
fluids of the same density, revolving as a rigid body. 

The case of Moon and Earth is well represented by a =^ 1, A = 4*333 ; for this 
gives 1 to 81 -35 as the ratio of the masses. With these values I find 

C=l- [2-46626] c~'- [2-45443] c-*- [2'56358] cr^ - [2-65073] c^^ - [2-72328] c"^, 

and 

r + R= 5-333 + [2-59898] c"^ + [2-24358] c"^ + [2-38748] c"^ + [2-77563] c^^ 

+ [3-32042] 0-7. 

Now c = 7-0 will be found to make C vanish, and, with this value of c, c — (r + R) 
= -414. 

If A be 4000 miles, c =: 6500 miles, and c — (r + ^) — 380 miles. 
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Thus, as far as this investigai.tion goes, it appears that when the fluid Moon is on the 
point of breaking up from stress of tidal and centrifugal forces the distance between 
the centres of Moon and Earth is 6500 miles, and the shortest distance between the 
two surfaces is 380 miles. 

This result must, however, from the native of the approximation, be an under- 
estimate of the distances. 

The whole of the present section has been suggested by a pamphlet by Mr. James 
Nolan "^ in which he criticises some of my previous papers. I have commented 
elsewhere on his criticisms.t 



§ 12. On the Case where the two Masses are unequal 

The results of the previous section point to a very remarkable limitation to the 
possibility of approach of two masses of unequal size. It has, therefore, seemed worth 
w^iile to consider this case numerically, and a case is therefore chosen which shall 
approach near to that which we know is the limit of possibility. I choose, therefore, 
a= 1, A = 3, which makes the ratio of the masses 1 to 27, and c = 5*3, which brings 
the protuberances into close proximity. 

The numerical details are omitted, but figs. 5 and 6 (Plate 23) give the results, the 
numerical values of the radius-vectors being, as before, entered on the figure. 

The elongation of the smaller mass is so- extreme that it is obvious that, rigorously 
speaking, the spherical harmonic approximation must be considered to break down. 
Nevertheless, I conceive that these curious figures may be held to indicate the general 
nature of the true result. 

It is remarkable that the smaller mass exhibits a marked furrowing round the 
middle. This seems to indicate that such a system tends to break up by the separa- 
tion of the smaller mass into two parts. 

§ IS. Summary, 

The intention of this paper is, first, to investigate the forms which two masses of 
fluid assume when they revolve in close proximity about one another, without relative 
motion of their parts ; and, secondly, to obtain a representation of the single form of 
equilibrium which must exist when the two masses approach so near to one another as 
just to coalesce into a single mass. 

When the two masses are far apart the solution of the problem is simply that of the 
equilibrium theory of the tides. Each mass may, as far as its action on the other is 

^ ' Darwin's Theory of the CTenesis of the Moon/ Robertson, Melbourne, Sydney, Adelaide, and 
Brisbane, 1885. 

t 'Kature,' February 18 and July 29, 1886. 
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concerned, be treated as spherical, and the tide-generating potential is given with 
sufficient accuracy by a single term of the second order of harmonics. As the masses 
are brought nearer to one another, this approximation ceases to be sufficient, terms of 
higher orders of harmonics become necessary to represent the potential adequately, 
and the departure from sphericity of each mass begins to exercise a sensible deforming 
influence on the other. 

When the departure from sphericity of one body produces a sensible deformation in 
the other, that deformation in its turn reacts on the first, and thus the actual figure 
assumed by either mass may be regarded as a deformation due to the primitive 
influence of the other mass, on which is superposed the sum of an infinite series of 
reflected deformations. 

But each mass is deformed, not only by the tidal action of the other, but also by its 
own rotation about an axis perpendicular to its orbit. The departure from sphericity 
of either body due to rotation also exercises an influence on the other, and thus there 
arises another infinite series of reflected deformations. It is shown in this paper how 
the summations of these two kinds of reflections are to be made by means of the 
solution of three sets of linear equations for the determination of three sets of 
coefficients. 

The first set of coefficients are augmenting factors„ by which the tides of each order 
of harmonics are to be raised above the value which they would have if the perturbing 
mass were spherical. It appears that, the higher the order of harmonics, the more do 
these factors exceed unity. 

The second set of coefficients correspond to one part of the rotational effects. They 
appertain to terms of exactly the same form as the tidal terms, and in the final result 
the terms to which they apply become fused with the tidal terms. These terms are 
the zonal harmonics of the several orders with respect to the axis joining the centres 
of the two masses. 

The third set of coefficients correspond to the remainder of the rotational effect, and 
they appertain to a different kind of deformation. These deformations are represented 
by sectorial harmonics involving cos 2<^, where ^ is azimuth measured from the plane 
passing tlirough the axis of rotation of the system and the centres of the two masses. 
That term of this set which is of the second order of harmonics, and which represents 
the ellipticity of either mass augmented by mutual influence, is the only term which 
is considerable, even when the two masses are very close together ; but the existence 
of the other harmonic deformations of this class is interesting. We may say, then, 
that all the tides of either mass are augmented above the values which they would 
have if the other mass were spherical ; that the ellipticity corresponding to rotation is 
augmented ; and that the deformation due to rotation is no longer exactly elliptic- 
spheroidal. 

The angular velocity of the system is found by the consideration that the repulsion 
due to centrifugal force between the tw^o masses shall exactly balance the resultant 

MDOCCLXXXVII. A. 3 H 
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attraction between them. If the masses were spherical, the square of the orbital 
angular velocity, multipUed by the cube of the distance between the centres, would be 
equal to the sum of the masses. When the masses are deformed, however, this law is no 
longer true, and the angular velocity has to be augmented by a factor a little greater 
than unity, which depends on the amounts of the deformations. 

The theory here sketched is applied above numerically to several cases, and the 
results will be found in the preceding paragraphs. We shall first consider the cases 
where the two masses are equal to one another. 

In the first example {J3 = y) solved numerically, the distance between the centres 
of the two masses is 2*83 times the mean radius of either of them. The two bodies 
are found to be elongated until they approach near to one another; but, as the 
character of the distortion is better illustrated in a subsequent case, the result is not 
given graphically. Ail the data, however, are found which will enable the reader to 
draw the figure if he should wish to do so. 

In the next example {fi = ^), with the masses still equal, the distance between the 
centres is reduced to 2*646 of the mean radius of either. The result of the solution is 
illustrated by two figures. In fig. 4, Plate 22, the section of the masses by a plane 
perpendicular to the axis of rotation is shown, and in fig. 5, Plate 23, we have the 
section by a plane passing through the axis and the centres of the two masses. On 
both figures are inscribed the values of the radii for each 15° of latitude in terms of 
the mean radius as unity, and the mean sphere, from which the distortion is computed, 
is marked by a short line on each radius. The elongation of the masses is, of course, 
considerably greater in the section through the axis than in the other section. Each 
mass is shaped somewhat like ati egg, and the small ends face one another and come 
very nearly into contact. 

In the headings to the figures^ amongst other numerical data, are given the square 
of the angular velocity and the angular momentum of the system. The density of 
the fluid being unity, the angular velocity co is given by the value of co^/iir ; this is 
the function of angular velocity which is usually given when reference is made to figures 
of equilibrium of rotating fluid, such as the revolutional or Jacobian ellipsoids of 
equilibrium. The moment of momentum of the system is given by reference to the 
angular velocity of a sphere, of the same mass as the sum of our two masses, rotating 
so as to have the same momentum. If, in fact, b be such a length that a sphere of 
fluid of that radius has the same mass as our system (so that b^ = ctP + ■^^)^ then the 
moment of momentum is given by a number (jl in the expression (f tt)' b^ X /x. By this 
notation the angular velocity and moment of momentum are made comparable with 
the results given in a previous paper'^^^ on the Jacobian ellipsoid of equilibrium. 
From that paper the following table of the axes, angular velocity, and moment of 
momentum of several solutions of Jacobi's problem is extracted. 

^^ ' Roy. Soo. Proc.', toI. 41, 1887, p. 319. 
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Jaoobi's Ellipsoids. 
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In figs. 4 and 5, Plates 22, 23, o)^J4.tt is '038, and the momentum /x is '472. On com- 
parison with the Table of Jacobi's ellipsoids^ we see that this corresponds with a con- 
siderably slower rotation than the 6th solution^ and nearly the same moment of 
momentum. 

In the next case the two masses are still closer (^ = \), the distance between the 
centres being only 2*449 times either mean radius. The result is illustrated in 
figs. '2 and 3 ; the explanation of figs. 4 and 5 serves^ mutatis mutandis, for these 
two figures also. 

This case is interesting because the masses have approached so near to one another 
that they partially overlap. Two portions of matter cannot, of course, occupy the 
same space, and the continuity of figures of equilibrium leads us to believe that the 
reality must consist of a single mass of fluid. In figs. 2 and 3 conjectural dotted 
lines are drawn to show how it is probable that the overlapping of the two masses is 
replaced by a neck of fluid joining them. The figures as thus amended serve to give 
a good representation of the single dumb-bell shaped figure of equilibrium. 

The angular velocity is here given by co^/Att = '049, and the moment of momentum 
by '482. In the sixth entry of the Table of Solutions of Jacobi's problem we find 
o)^/Att = '0536, and the moment of momentum /x = '481. This ellipsoid has, then, the 
same moment of momentum, and only about 4 per cent, more angular velocity, than 
our dumb-bell. It has seemed, therefore, worth while to mark (in chain-dot) on 
figs. 2 and 3 the outline of this Jacobian ellipsoid of the same mass as the dumb- 
bell. The actual vertex of the ellipsoid just falls outside the limits to which it was 
possible to extend the figure. 

In the paper above referred to it is shown how the enei-gy of the Jacobian ellipsoid 
is to be computed. If we denote the kinetic energy by (Itt)^ b^ X e, and the intrinsic 
energy by (f7r)^b^ X [i — 1),* then it appears that in the case of the ellipsoid drawn 
in these figures e = -0964, i = '4808, and tlie total energy B =^ e -{' i := '5772. 

* The intrinsic energy being negative, it is more convenient to tabulate i a positive quantity. 

3 H 2 
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Now in the case of our dumb-bell figure it appears, from calculations referred to in 
the Appendix, that e~ '0925, i= '4873, and E= '5798. Hence in the dumb-bell 
figure the kinetic energy is less, but the intrinsic energy is so much greater that the 
total energy is about a half per cent, greater. These numbers are^ of course, computed 
from the approximate formulae, and must not be taken as rigorously correct for the 
dumb-bell figure of equilibrium. 

With reference to a figure of transition from the Jacobian ellipsoid. Sir William 
Thomson has remarked : — ^'^ 

'^ We have a most interesting gap between the unstable Jacobian ellipsoid, when 
too slender for stability, and the case of smallest moment of momentum consistent 

« 

with stability in two equal detached portions. The consideration of how to fill up 
this gap with intermediate figures is a most attractive question, towards answering 
which we at present oflPer no contribution." t 

Figs. 2 and 3 are intended to form such a contribution, but it is certain that the 
matter is far from being probed to the bottom. 

M. PoiNCARE has made an admirable investigation of the forms of equilibrium of a 
single rotating mass of fluid, and has especially considered the stability of Jacobi's 
ellipsoid. J He has shown, by a difficult analytical process, that when the ellipsoid is 
moderately elongated (he has not arrived at a numerical result) instability sets in by a 
furrowing of the ellipsoid along a line which lies in a plane perpendicular to the 
longest axis. It is, however, extremely remarkable that the furrow is not symmetrical 
with respect to the two ends, and thus there appears to be a tendency to form a 
dumb-bell with unequal bulbs. 

If M. Poincarje's result shall appear to be not only true, but to contain the whole 
truth concerning the mode in which instability of the ellipsoid supervenes, then there 
must be some other transitional form between the unsymmetrically furrowed Jacobian 

"^ Thomson and Tait's 'Natural Philosopliy,' 1883, § 778" (i). 

f In 778" (g) he remarks that ''a deyiation from the ellipsoidal figure in the way of thinning it in 
tlie middle and thickening it towards the end would, with the same moment of momentnm, give less 
energy." I conceive that the energy referred to throngliout this paragraph is kinetic only, and we have 
seen that tke kinetic energy is less for the dnmb-bell than for the ellipsoid. 

[If we write U for a quantity proportional to the excess of kinetic above intrinsic energy, so that 
U = 6 + (1 — ^), then fiigures of equilibrium are to be determined by making U stationary for variations 
of the parameters involved in it. This course is actually pursued in tbe Appendix below, the function 
(viii.) being, in fact, this U ; and the variations of it, being made stationary, afford a controlling solution 
of the problem of this paper. The similar method may easily be applied to the ca.se of Jacobi's ellip- 
soids. From this point of view the interesting function to tabulate is e -f (1 — -i^), and we observe that 
in the case of the Jacobian ellipsoid referred to on the last page it is '6052, and for the dumb-bell it is 
•6156. Is not tbe energy referred to by Sir W. Thomson this function U ? (Addition to foot-note^ 
dated October 10, 1887.)] 

J ' Acta Mathematica,' 7, 3 and A, 1885. 
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ellipsoid and the dumb-bell; except, perhaps, in the case where the two bulbs pass on 
to two masses of a definite ratio. 

M. PoincaFvE's work seemed so important that this paper was kept back for a year, 
whilst I endeavoured to apply the principles, which he has pointed out, to the dis- 
cussion of the stability of the two masses. The attempt, which is given in the 
Appendix, is apparently abortive, on account of the imperfections of spherical harmonic 
analysis when applied to bodies which depart considerably from the spherical shape. 

We must, therefore^ leave this complex question in abeyance^, and merely point to 
the Appendix as an example of the method which must almost certainly be pursued 
if this problem is to yield its answer to analysis. 

Allusion has just been made to the imperfection of spherical harmonic analysis, and 
this brings us naturally to face the question whether that analysis may not have 
been pushed altogether too far in the computation of the figures of equilibrium under 
discussion. This question is considered in § 9, and a rough criterion of the limits of 
applicability of this analysis is there found. From this it appears that even in the 
cases of figs. 2 and 3 the result must present a fair approximation to correctness. 
The criterion, indeed, appears to be such as necessarily to give too unfavourable 
a view of the correctness of the result. 

The rigorous method of discussing the stability of the system having failed, certain 
considerations are adduced in § 11 w^hich bear on the conditions under which there is 
a form of equilibrium consisting of two fluid masses in close promixity. It appears 
that there cannot be such a form with the two masses just in contact, unless the 
smaller of the tw^o masses exceeds in mass about one-thirtieth of the larger. 

If we take into consideration the fact that the criterion of the applicability of 
harmonic analysis is too severe, it appears to be worth while to find to what results 
the analysis leads when two masses, one 27 times as great as the other, are brought 
close too'ether. The numerical work of the calculation is omitted, since the numbers 
can only represent the true conclusion very roughly ; but the result is illustrated 
graphically in figs. 6 and 7, Plate 23. These figures can only serve to give a general 
idea of the truth, but the form into which the smaller mass is thrown is so remarkable 
as to be worthy of attention. The deep furrow round the smaller mass, lying in 
a plane parallel to the axi^ of rotatiqn, cannot be due merely to the imperfection 
of the solution ; and it appears to point to the conclusion that there is a tendency for 
the smaller body to separate into two, just as we have seen the Jacobian ellipsoid 
become dumb-bell shaped and separate iixto two parts. 

In this paper, indeed, we have sought to trace the prqce^s in the opposite direction, 
and to observe the coalescence of two masses into one. The investigation is comple- 
mentary to, but far less perfect than, that of M. Poincar^, who describes the series 
of changes which he has been tracing in the following words : — 

'' Considerons une masse fluide homogene animee originairement d un mouvement 
de rotation ; imaginons que cette masse se contracte en se refroidissant lentement, 
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mais de fagon a rester toujours homogene. Supposons que le refroidissement soit 
assez lent et le frottement interieiir du flnide assez fort pour que le mouvement de 
rotation reste le m^me dans les diverses portions du fluide. Dans ces conditions le 
fluide tendra toujours a prendre une figure d'equilibre seculairement stable. Le 
moment de la quantite de mouvement restera d'ailleurs constant. 

'' Au debut, la densite etant tres faible, la figure de la masse est un ellipsoide de 
revolution tres peu different d'une sphere. Le refroidissement aura d'abord pour effet 
d'augmenter I'nplatissement de rellipso'ide, qui restera cependant de revolution. 
Quand Taplatissement sera devenu a peu pres egal a f , Tellipsoide cessera d'etre de 
revolution et deviendra un ellipso'ide de Jacobi. Le refroidissement continuant, la 
masse cessera d'etre ellipsoi'dale ; elle deviendra dissymetrique par rapport au plan 
des yz, et elle aflfectera la forme representee dans la figure, p. 347.^^ 

^^ Comme nous Tavons fait observer a propos de cette figure, I'ellipsoide semble se 
creuser legerement dans sa partie moyenne, mais plus pres de Tun des deux sommets 
du grand axe ; la plus grande partie de la matiere tend a se rapprocher de la forme 
spherique, pendant que la plus petite partie sort de Tellipsoide par un des sommets 
du grand axe, comme si elle cherchait a se detacher de la masse principale. 

'* II est difficile d'annoncer avec certitude ce qui arrivera ensuite si le refroidisse- 
ment continue, mais il est permis de sup[)oser que la masse ira en se creusant de plus 
en plus, puis en s etranglant dans la partie moyenne, et finira par se partager en deux 
corps isoles. 

'' On pourrait etre tente de chercher dans ces considerations une confirmation ou une 
refutation de I'hypothese de Laplace, mais on ne doit pas oublier que les conditions 
sent ici tres differentes, car notre masse est homogene, tandis cj[ue la nebuleuse de 
Laplace devait etre tres fortement condensee vers le centre.'' t 

It was in the hope that the investigation might throw some light on the nebular 
hypothesis of Laplace and Kant that I first undertook the work. It must be 
admitted, however, that we do not obtain much help from the results. It is justly 
remarked by M. Poincare that the conditions for the separation of a satellite from a 
nebula differ* from those of his problem in the great concentration of density in the 
central body. But both his investigation and the considerations adduced here seem 
to show that, when a portion of the central body becomes detached through increasing 
angular velocity, the portion should bear a far larger ratio to the remainder than is 
observed in the satellites of our system as compared with their planets ; and it is hardly 
probable that the heterogeneity of the central body can make so great a difference in 
the result as would be necessary if we are to make an application of these ideas. It 
appears then at present necessary to suppose that after the birth of a satellite, if it 
takes place at all in this way, a series of changes occur which are quite unknown. 

* The furrowed ellipsoid of Jncort, 

t Poincare, ' Acta Mathemat.,' 7, 1885, p. 379. 
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Appendix. 



On the Energy and Stability of the System. 



M. PoiNCARE bas shown in his admirable memoir, referred to in the Summary, how the dynamical 
stability of a rotating fluid system in relative equilibrium dejDends on the energy. Certain factors in 
the expression for the energy, which he calls coefficients of stability, are there proved to afford the 
required criterion. 

It will now be shown how in this case these coefficients of stability are determinable, at least as far as 
spherical harmonic analysis permits* The results will also cast an interesting light on the methods by 
which the equations to the two masses have been obtained. 

The task before us is to determine the " exhaustion of potential energy " of the two masses in 
presence of one another as due to the deformation of each from the spherical figure by yielding to 
gravitation and to centrifugal force. 

The work will be rendered simpler by the introduction of a new notation. Let us write, then, as the 
equations to two shapes, which are not necessarily together a figure of equilibrium : — 



a 

E 
A 



1 4. ^'1'' 2/^+ I /J\3 /a\^-3 
jc=2 2k — 2\o / \c 

. %'^2k + l/a\^ fAY'^ 
■^ 1Z2 2k ^ 2 



nk 



lOjc 



r 



k 



:Pk 



^^Wk + o. 



/yt. 



k 



1 ^ E' 



P 



k 



^^TT; 



> 



>^ + 2 



i# 



-J 



(i.) 



It will be observed that these equations have the same form as (72), but that the constants introduced 
are different from the h, I, m, e, which were determined, so that the figures might be figures of equilibrium. 
At present we do not assume that (i.) do represent figures of equilibrium. 

The energy lost may be divided into several parts : — 

e^, the energy lost by the mass a yielding from its spherical figure to the gravitation of the mean 
sphere a. 

eg, the exhaustion of mutual energy of that layer of matter on the mass a which constitutes its 
departure from sphericity. 

63, the loss of energy due to the deformation of the mass a in presence of the mean sphere A. 

JE-^, E^, ^3, the similar quantities for the mass A. 

(JEe)^, the loss of mutual energy of the two layers in presence of one another. 

65, the loss of energy due to the deformation of the mass a in the presence of centrifugal force duo 
to rotation w. 

JEJ^, the similar loss for A, 

1st. Gi is equal and opposite to the work required to raise each element of the layer on a through half 
its own height against the gravity due to the mean sphere a. This gravity is ^Tra. Co-latitude and 
longitude being denoted by 0, 0, let dixr = sin dO dcj), an. element of solid angle. In effecting the inte- 
grations, the properties of spherical harmonic functions are used without comment, viz. : — • 



W/c 



r 



k 



dw = 



^TT 



2^ + i ' 



'l^^Wk- 



3 



.^k 



dw = 



47r k-\-2\ 



I 



ft/ <j 



2k ^ I 2.k~'2\ 



lUk^^Wk+Qdw = 0, 



WiWk dw = 0, 



Wih-wjc..odw = 0. 



Then, taking only a typical term of the first of (i.) 

r5 1 



i-7ra' 



• 2 



2k + 1 MV /a^^-2 



2k -2 \c 



n 






'IVk 



,^2, 



-^^k ^^ . ^k 



■ /yK I 



"1 Q 



dw 



iTT 



3 



\o^ 



{2k ^2y \c) [c^ 



ni^ + 



k-\-2\ 
2,k-~-2\ 



19k' 
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wbence, with all the terms, and i^ememberiiig that (a/c)^ = <y, (iL/c)^ = F, 



e-L = 



1 

2 






c 



_3 

2 



'^V^ 



= CO n 



'^^ + 1 



Jc-l r 



0/ /.■ = 2 



,.==2 2/^-2 /^ 



2 . /c — 2 



(ii.) 



The formula for TJi m.ay be written down bj symmtetrj. 

2nd. ^2, the exhaustion of mutual energy of the lajer on itself, is half the potential of the layer afc 
any element, multiplied by the mass of the element, and integrated over the w^hole sphere. 

The potential of the layer is 



W3.J» 3 



a\^ /a^^+^ 



3c ^- = ^^i ~~ 2\ol \7' 






^hvk+2 



Then, at an element of the layer r = a, and taking a typical term only, we haye 



^2 = 



_i 4.7rA^ 3(2Z^ + l)/aY^^~V^V^ 3 



o. • 



3o ' (2/^-2)3 \oJ \c 



^^,(^^J + p,(^ 
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Jc 



dw^ 



whence 
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■47r\3 jBa^' 
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3 



1 



k-i 



0/ k = 2 



: = o2k-2 h 



njc^ + 



fe + 2 1 

2T¥^2! 



pr 



(iii.) 



The formula for JEJg may be written down by symmetry. 

The addition of e^ to e^, and of JS^ to ^3, simplifies these expressions by cutting out the factor 
immediately following the X iiJ. either, and replacing it by unity. 

3rd. 63 is the loss of energy due to raising the layer on a in presence of the mean sphere A, Wo 
multiply the potential of the sphere A by the mass of the element on a, and integrate throughout angular 
space. 

The potential of the sphere -4, when transferred to the origin 0, is 

4i7rA^ ^'T^'^ la^ It\^ ivjc 



3c lc = \G 



a I T 



h 



Then, at an element of the layer r = a and taking a typical term, 



47r43 2/^^ + 1 lA\^la^^-^ 



- — 



3c 2/i; ■— 2 \ c / \ c 



w- 



Im' 



*} J 



"*(5) '^'^' 



whence 






njc. 



• * • • 



• • (iv.) 



The expression for i?g may be written down by symmetry. 
and (iv.), we have 



On collecting results from (ii.), (iii.), 



1^ ^u9\ 
^nr — —- :r; pr K ... (v.) 



4. ^--21 



and a similar expression for 'E^ -f- ^3 + ^/3. 

4th. (-E/e)^^ is the loss of energy of one layer in the presence of the other. We take the potential of 
the layer on A^ multiply it by the mass of an element on a^ and integrate. 

The potential of the layer on A when transferred to the origin 0, as in (22-ii.), is 



. 3 



3c ^ 



2 2 



.c 



'h -{- i\ V^~^ f(i\^ (t\^ ivjc 



h -\- i\ r*~^ _^ (o\^^ /r\^ ^hu/c^, 



^'-2 i = 2 L ^'1 '2^1 ^--1 \Gj \al r^ ^ — 21 /i? -h 2! f~ 1 



-P. 



.c \a 



+ 2 



Jc 



Introducing this into the integral, only taking a typical term, and neglecting those terms in the integral 
which must vanish, w^e get 
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air > . 



^, 



k 



Then, effecting the integrations, and putting in the S2, we get 
(-E?e)4 = 



}• 



V,f^ \. (vi.) 



This involves the two figures symmetrically. 

5th. 65 is the loss of energy in the yielding of the figure a to centrifugal force. To find it, we 
multiply the potential due to rotation by the mass of each element of the layer a, and integrate. 

By (35) and (52) we know that the rotation potential is 



^-w^ct^ 



r\3 rm 



Ch> 



f ^3 _ X ^H 1 



As this only involves harmonics of the second order, we may neglect in the layer a all terms except 
those of the second order. Thus we get 
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-n ,s 



} 



aw 



= 4*7r , -:^^w'^a^ 



4! 



4<wYA^oP f 3&^3/^\2 






IOttIc, 



K + 6p,) } . 



»»f»»..»«» I VXl. i 



The expression for .Ug may be written down by symmetry. Collecting results from (v.), (vi.), 
and (vii.), we get, for the whole exhaustion of energy, 
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• • 



(viii.) 



The expression is found without any assumption that the two masses are bounded by leYel surfaces, 
and therefore ip equilibrium. But the condition for equilibrium is that the difierential coefficients of 
B with respect to any one and all of the parameters % j?, JV, P, shall vanish. If we equate to zero 
dVldnk, we get 



1 —fik + 



^^^3i:^^00 p-1 J^^il 



2 i — 1 hlil 



Ni=0, 



If, howpver, ^ = 2, there is on the left-hand side an additional term 



Sw^ /a\^ ./AYn w^/c\^ 



3/tll 1- 



I6r \o/ ■' ^ \c/ 1 87r\A 

The equation of dV/dNi to zero gives a similar equation. 
If we equate to zero d Vld;pjc, we get 
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If, howeyer, ^ = 2, there is on the left-hand side an additional term 



Utt^VcI ' Mc/12!0! SttI^/ 4.3 "'^U 



The equation of dVjdPi to zero gives a similar equation. 

How, if we pnt hk + ^;5: for n/c, except when ^ = 2, and then put h^ + h + t\^ (cjAy = %, and 
similarly introduce the H'b and L*s ; and if we put pk = '^/^j except when ^ = 2, and then put 
p^ = '^^3 + yV^ (<3/^)^3 and similarly introduce the Jf' s, it is easy to see that the equations (i.) to the 
two surfaces become the same as (72), and the equations of condition between n and N, and between 
p and P, become exactly those which w^e found by a different method above in (23), (44), and (59). 
The only difference is that the equations for h and I are fused together. 

This, therefore, forms a valuable confirmation of the correctness of the long analysis employed for the 
determination of the forms of equilibrium. 

The formula (viii.) also enables us to obtain the intrmsic energy of the system, that is to say, the 
exhaustion of energy of the concentration of the matter from a state of infinite dispersion to its actual 
shape, with its sign changed. 

The last line of (viii.) depends on the yielding of the fluid to ceDtrifugal force, and must be omitted 
from the exhaustion of energy. 

Besides the rest of (viii.), we have in the exhaustion of energy of the system, the exhaustion of the 
two spheres and their mutual exhaustion. 

It is clear, then, that the intrinsic energy is 







c I h -— 1 






k=2 i=2 

where the w, N, P, p, have their values determined in accordance with the condition that the surfaces are 
level surfaces. 

In evaluating the intrinsic energy from this formula, it is convenient to refer the energy to that of a 
sphere of such radius, b, that its mass is equal to the whole mass of the system. Then h^ =z a^ -]- A% 
and we may take the intrinsic energy as 

(|^)H5(*-l). 

Thus i will be a numerical quantity which is positive. 

I find from (ix.), with /5 = ^ and a = A, that t = '4873. 

With regard to the kinetic energy of the system, we have seen in § 10 that the moment of momentum 
is (|7r)t6^X/4, where /t is a numerical quantity, and we find in the course of the determination the 
function 3w2/47r. Then, since the energy is half the moment of momentum multiplied by the angular 
velocity, it is clear that the kinetic energy is 



(1^)3 &5x|^ 



47r 



The kinetic energy, as represented by e = ^/Ji (3w2/45r)t, is comparable with the intrinsic energy as 
represented by 4. 

In the case of 13 =: i, and a — A,l find the kinetic energy e = '0925. 
Thns the total energy B = 6 + ?* = *5798.* 



* 



[See foot-note added to the Summary above, October 10, 1887.] 



EQUILIBRIUM OF ROTATUSTa MASSES OF FLUID. 



427 



If the energy of a system be expressed as the sum of a nnmber of coefficients, each multiplied by the 
square of a parameter, it has been shown by M. Poincaee that the stability of the system depends on the 
signs of these coefficients, which he calls " the coefficients of stability." But, if the expression for the 
energy involves the products as well as the squares of the parameters, the coefficients of stability are the 
roots of a determinantal equation involving the second differentials of the energy with respect to the 
parameters. 

Let F, a linear quadratic function of x, y, z, &c., be the energy of a system in equilibrium ; then the 
determinantal equation is 



dW 
dy dx' 

dW 
dz dx 



X, 



dW 
dxdif 

dW 



dy'^ 

dW 
dz dy 



-X, 



d?V_ 
dx dz 

dW 
dy dz 



&c. 



&c. 



dW 



dz^ 



&c., 



&c., 



' — X, &c. 



&o. 



= 0. 



(x.) 



The solution of this equation in X involves the determination of the several fundamental modes of 
vibration of tbe system ; and the roots are the coefficients of stability. 

Now suppose that F involves a constant, then, in causing that constant to vary continuously, we have a 
series of systems of equilibrium of the same kind ; and the coefficients of stability vary continuously at 
the same time. If the system be initially in stable equilibrium, the stability ceases when a coefficient of 
equilibrium vanishes. The system at the moment of instability is in a condition of " bifurcation," that 
is to say, there is another series of shapes of a different kind, of which this shape is a member. In 
making the constant vary past the critical value, we find this second series of shapes stable, whilst the 
first is unstable. 

If the system be in uniform rotation, so that instead of absolute equilibrium there is equilibrium 
relatively to uniformly rotating axes, the same theorems hold true, provided that only one root of the 
determinantal equation vanishes at a time. 

This last is the case which we are considering, and the constant, which we suppose to vary continuously, 
is c, the distance between the two centres of the mean spheres of radii a and A. 

When the two masses are far apart the equilibrium is stable, but when they are brought closer a time 
may come when one of the coefficients of stability vanishes. 

The condition for the vanishing of a coefficient of stability is determined by the determinant (x.) with 
X = 0. 

To find the determinant, we have to evaluate the second differentials of B with respect to the 
parameters w, JV, p, P. 

If we form the determinant corresponding to (x.) with X = 0, it is obvious that two infinite squares 
of entries which are diagonally opposite to one another, and which meet at a corner, are to be filled in 
with differentials involving dndp, dndP, dNdP, dN dp, in the denominators. All these entries are 
zero, and hence the infinite determinant splits into two independent infinite determinants, one only 
involving the differentials with respect to N, n, and the other only those with respect to P, p. The JV, n 
determinant may be called "the tidal determinant, '* fche P, p one "the rotational determinant"; for the 
origin of the terms in each is obvious. 

By considering only the tidal determinant, we see how the other may be treated very shortly. 

For the sake of brevity, write 

d^EldnjcdNi 



nk Ni = .' 



{--d^Eldnk^y^ (- d^Eldm)l ' 



(xi.) 



Then stability vanishes, as far as regards the tidal forces, when 

3 I 2 
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(xii.) 



It will be observed that, by tlie notation (xi.), and tlie appropriate division of tlie columns and rows of 
tlie tidal determinant, it has been converted into a determinant in whioli the diagonal consists of ones. 
If we drop certain factors common to the whole, we have, by differentiating (viii.), 



dn]? 
d^B 



3_ 
2 



^\8 ^Ic 1 



7 



G 



h-l 



dm 



2 



2 



'aV r^-i 



d^E 



d^B 



h-l 



dfik dni dN/c dMi 



0, 



dfik dNi 
dW 



= (l) 



^lA'fla\^ 7^-1 r^-i h + i\ 



.c 



k — 1 i 



dpk^ 

d^B 
dpkdFi 



3 
2 



ay 



c) 2.k~-2l h 
'A\^lo\^ ^ + 21 



d^B 



1 dFj? 

,1c — \ pi — 1 



a\3 it + 2! r^-i 



d^B 



d^B 



r 



2.h-2l h 






c/ 2.h- 
h -\- i\ 

l + il 



2! h 



d;pjc djpi dPjc dP 







2 ! /b - 2 ! ' 



(xiii.) 
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With these values (xiii.) we easily find, by substitution in (xi.), that 
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[This expression gives the value of each of the entries in the infinite determinant (xii.). 

E'ow it is possible, by a certain laborious investigation which I do not here reproduce, to develope 
this infinite '* tidal " determinant in the form of a series, and then to show that, however close the two 
masses may be to one another, the series arising from the tidal determinant can never vanish. It may 
also be proved that the other infinite determinant, which results from "rotational" terms, is necessarily 
greater than the tidal determinant, and a fortiori can never vanish.]* 

Thus it might be held that stability must subsist in the figures of equilibrium until the two masses 
come into contact. But, as appears from § 11, it is certain that, if one of the masses be smaller than the 
other, this cannot be the case. In fact, the investigation must break down on account of the imperfection 
arising from the use of spherical harmonic analysis. 

We have seen that the infinite determinant, which gives the coefficients of stability, splits into two 
parts when we rely on spherical harmonic analysis ; but when instability ensues it must be brought 
about by the joint action of the tidal and rotational forces. It appears certain then that, if a rigorous 
analysis were used, this separation would not take place. 

Although the present investigation proves thus to be abortive, I have thought it worth while to sketch 
the process, because it almost certainly indicates the line that must be pursued whenever a more 
rigorous analysis shall be applied to this difficult problem. 



* This paragraph, replacing the investigation referred to, added July 12, 1887. 
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